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SUMMARY 


The actual buckling stress cr cr can be calculated from the first 
author's exact theory as well as by his method of split rigidities. Both 
methods yield practically identical results. By the latter method simple 
formulas are obtained which express the actual buckling stress a cr 
directly in terms of the column and local or plate buckling stresses. 
Columns with box, I-, H-, and T-sections and angles are considered 
separately. Interaction of practically significant magnitude occurs o nl y 
in cases of flexural and torsional buckling. In these cases the addi- 
tional effect of distortion of the cross section is also taken into 
account. The theory includes buckling in the plastic range. No post- 
buckling phenomena are considered in the theoretical part of the paper. 

Tests were carried out for a considerable range of ratios of cor- 
rected free length to radius of gyration on two sections, for one of 
which the local buckling stress was in the plastic domain, and for the 
other, in the elastic domain. The experimental buckling stresses are 
in excellent agreement with those predicted by the theory. 


INTRODUCTION 


It is customary to consider that a column may buckle in either one 
of two ways: (a) By deflection of the entire column in a half wave of 

length equal to the effective column length (column buckling) or (b) by 
plate buckling of its component webs and flanges in shorter or longer 
half waves (local or plate buckling). In the first case it is tacitly 
assumed that no distortion of cross section occurs, while in the second 
the lines of intersection of the midplanes of the various plates are 
assumed to remain straight. For a given column, buckling is supposed 
to occur at the lower of the two critical stresses, column or local. 

In reality, however, there is an interaction of these two modes of 
buckling, so that the real buckling stress cr cr will be smaller than 
either of the buckling stresses for column or local buckling. 
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With column buckling the buckling stress Oj is determined by the 
Euler or Engesser load. In figure 1 0 -^ is plotted against the ratio 

P = a/b of the half wave length a to the web width b, for example, 
for a column section like that in figure 2(a). On the other hand, with 
local buckling, which assumes the lines of intersection of the middle 
planes of the plates to remain straight, the buckling stress is given 
by a 2 in figure 1. The latter becomes minimum for a ratio P x = a/b 
of order of magnitude 1. If no interaction is taken into account for 
p < P 2 in figure 1 the minimum plate buckling stress ( g p) TT n , n is 

smaller than o^. Hence for an I-section with an effective buckling 
length a, where a/b < P g , the plate buckling stress a 2 is governing 
and web and flanges will buckle in relatively short waves. If a/b > P 2 
the column buckling stress cr-^ governs and the column buckles as a 
whole in a single half wave. 


Actually the buckling deflection consists of a deflection w^ of 
the member as a whole, as it occurs with column buckling, and a deflec- 
tion w 2 of the web, as it occurs with plate buckling (fig. 2(b)). ^ 
Assuming first an infinite rigidity against column buckling, w-l will 
be zero and a buckling of the web will occur at a stress a 2 with a 
maximum deflection w 2 and in waves with a half wave length a of the 
order of magnitude of the web width b (a/b = Pp in fig. l). This web 
deflection w 2 will cause an external moment in the column as a whole. 

It will result in an entirely negligible internal moment in the column, 
however, since the latter is practically exclusively caused by the 
deflection of the flanges alone. Hence, assuming the column again to 
have a finite rigidity, this deflection w 2 has the same effect as an 
initial deflection w^ < w 2 of the entire column, causing an external 
but not an internal moment. It is well-known that an initial deflec- 
tion w i of a column causes an extra deflection of about |?cA °E “ a C'J w i> 
where o c is the actual compressive force and Cg is the elastic buckling 
stress of the column (see reference l) . Since, for the very small half 
wave length a of the deflection w 2 , the column buckling stress a ± 


for multiple-wave buckling is in most cases very high as compared with 
the plate buckling stress a 2 (compare ax and a 2 for a/b = Pp in 
fig. 1), this means that the deflection w 2 causes a column deflection w x 

which is smaller than [a 2 /(ap " a S\ w 2 and hence very sma11 aS compared 
with w 2 » This extra deflection w-^ increases the deflecting forces 

-ta dx dy (d 2 w/dx 2 ) that act on an element t dx dy of the web, while 

it increases only slightly the restraining forces acting on that element. 
Since w^ « the increase of the deflecting forces will be relatively 


NACA TN 2640 


3 


small, so that the actual buckling stress a QT will be only slightly 
less than a 2 (compare a cr and cr 2 for a/b = 0-]_ in fig. l). 

If, on the other hand, the I-section has an effective buckling 
length a which is larger than the limiting value a 2 = p 2 b, so that 

P = a/b is larger than p 2 in figure 1, buckling will occur in a single 
half wave. In this case, besides the deflection w-^ of the member as a 
whole, a small deflection w 2 of the web (fig. 2(b)) will also occur 


element t dx dy. 

In this case, in a similar way as above, the extra deflection w 2 
will be of the order of magnitude [hi/(a 2 - a-jJ w^_ which is small as 

compared with w^ because for a/b ^ P 2 , in most cases, cr 2 is very 
much larger than (fig. l) . The deflection w 2 increases the 

external moment in the column while it practically does not increase its 
internal moment. Since, however, w 2 is small with respect to w^, the 
actual buckling stress a cr will be only slightly smaller than a^. 

The qualitative conclusions above are fully worked out quantitatively in 
the theoretical part of the paper and are confirmed, in particular, by 
study of the interaction equations (74) and (75). 

Considerable interaction occurs and governs the actual buckling 
stress only if the individual buckling stresses and a 2 have their 

smallest value for a half wave length about equal to the effective column 
length. This happens, for example, with angles and T-sections, where for 
flexural buckling as well as for torsional buckling (which is a simplified 
form of plate buckling) both individual buckling stresses and cr 2 

are the smaller the larger the length of the column. In the latter case 
interaction actually occurs between three individual modes of buckling 
(fig. 3). 

In the theoretical part of the paper the exact theory for determining 
the actual buckling stress cxcr of arbitrary plate assemblies is given. 

It leads to transcendental buckling conditions for cr cr which have to be 
solved by trial and error. Furthermore, by the first author's method of 
split rigidities design formulas are derived which express a cr directly 
in the individual buckling stresses cr^ and o 2 and for a column with 
T-section in a 1} a 2 , and a^. As is customary, in all derivations 

column deflections by shear were neglected. 


In the experimental part of the investigation stress-strain tests 
were carried out on square tube columns. Column tests were then carried 
out using the stress-strain data to study the interaction of local and 

C Cl I nmvi "W, 1 - "1 -J 


■because of the deflecting forces 



that act on each 


column buckling. 
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This investigation was carried out at Cornell University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. The theoretical part was carried out by 
Professor P. P. Bijlaard and the experimental part was carried out by 
Professor G. P. Fisher. The project was directed by Professor George 
Winter, who is also at Cornell University. 


SYMBOLS 


a 


L 

b 

c,k 

e = (E/E s ) - 1 


half wave length of buckling 
column and plate length 
plate width 
constants 


n = s + 6qr 
P 

p y 

q = a-,2 _ (bA. 2 /d) 

q = 0^2 _ VX. 2 

r = a 2 2 + (b\ 2 /d) 

r = a 2 2 + VX, 2 


number of waves (used in section entitled "Exact 
Theory") 


equivalent load, including influence of twisting 
moments, transferred by plate to beams 

in plastic range (used in section entitled "Exact 
Theory" ) 

in elastic range (used in section entitled "Exact 
Theory") 

in plastic range (used in section entitled "Exact 
Theory" ) 


in elastic 
Theory" ) 


range (used in section entitled "Exact 


radius of gyration 


polar radius of gyration about shear center 


O 

3 1,2 = ]Fq) (®1, 2^ ~ A l,2 a cr) in plastic range 
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%2 * T 1 ?^ 2 ■ 

- A]_j2 a cr) in elastic range 

t 

plate thickness 

u = S - 0r 2 

(used in section entitled n Exact Theory") 

v = S - 6 q 2 

(used in section entitled "Exact Theory") 

w 

deflection 

x,y,z 

coordinates 

y 0 

distance between shear center and center of gravity 

A,B,D,F,G,H,K 

constants in theory of plastic plate buckling (used 
in section entitled "Exact Theory") 

A l,2 

cross section of beams (used in section entitled 
"Exact Theory") 


total cross section of columns 

K 

total cross section of webs situated perpendicular 
to direction of buckling 

Af 

total cross section of flanges situated perpendicular 
to direction of buckling 

B l,2 

flexural rigidities of beams (used in section entitled 
"Exact Theory") 

c 1 ,c 2 ,c 3 ,c 1| . 

constants 

C w 

warping constant 

D 

deflecting force in method of split rigidities 

E 

modulus of elasticity 

E s 

secant modulus 

E t 

tangent modulus 

G 

modulus of rigidity 

I 

moment of inertia 


6 


NACA TN 26 i+0 


M 

N 

Q 

R 

S 


polar moment of inertia about shear center 
Saint Venant torsion constant 
bending or torsional moment 
plate flexural rigidity 
transverse shear 

resisting force in method of split rigidities 
shear center 


V 

X,Y,Z 


a 


l,2> a l,2 
3 

°-)7 y £ 

n 

X = pjf/L 

V 

CT 

t = q + r 

+ = V tCT cr/ N 


energy 

coordinate axes 

given by equations (8), (ll), ( 27 ), and (28) 
ratio between half wave length a and plate width b 
coefficients in method of split rigidities 
reduction coefficient for plasticity 
(used in section entitled "Exact Theory") 

Poisson's ratio 
normal stress 

(used in section entitled "Exact Theory") 

(used in section under "Buckling of I -Section in 
Direction Perpendicular to Plane of Web" entitled 
"Comparison with Exact Theory") 


i = Ip W /I p , where I p is polar moment of inertia of web about shear 
center (used in section under "Buckling of Columns with 
T-Section" entitled "Critical Stress for T-Section with 
Fixed Shear Center Axis") 

^ reciprocal spring constant of restraining plate 

(used in section entitled "Exact Theory") 
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G = EIDi if 
0 = /ta cr /EI 


(used in section entitled f, Exact Theory”) 
(used in section entitled "Exact Theory”) 


THEORETICAL INVESTIGATION 

EXACT THEORY 


The exact general theory of elastic as well as plastic buckling of 
plate assemblies was published by the first author (references 2 and 3). 
It will be summarized here to the extent that it is used in this paper. 

Consider an asymmetric column, consisting of a web plate with 
width b and thickness t, which is supported at both unloaded edges 
by flanges of different width 2b ! and thickness t T (fig. 4(a)). Let 
the web plate be the "buckling” plate which is restrained at both edges 
by the flanges, the "restraining" plates. 

The X-axis is chosen in the longitudinal direction of the column, 
while for the buckling plate the Y- and Z-directions are chosen as shown 
in figure 4. Hence the differential equation for the buckling plate is, 
for the general case of plastic buckling (references 3 and 4), 


El 


A ^ + 2 (B + 


2F V 


A 


+ t o 


cr 


dx^ 


= 0 


( 1 ) 


where I = t3/l2 and a cr is the buckling stress cr x . For pure com- 
pression, as occurs here, A, B, D, and F, if expressed in terms of 
the secant modulus E s , the tangent modulus E^, and Poisson 1 s ratio v, 

(reference 4): 


are 
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A = Vi/% 

B = \|r 2 /% 

D = 

tl = 1 + 3(E-t/E s ) 

♦ 2 = 2 - 2(1 - 2V)(E t /E) 

t 3 = 4 

\ = (5 - 4V + 3e) - (1 - 2v) 2 (E t /E) 
F = l/(2 + 2 v + 3e) 
e = (e/E s ) - 1 


( 2 ) 


With v = 0.5 these values reduce to those used by Stowell (reference 5) 


A = (l/3)(E s /E) + (Et/E) 
B = (2/3)(E s /E) 

D = (4/3)(e s /e) 

F = (1/3)(E S /E) 


In the elastic domain in equations (2) e = 0 and E s = E^ = E, so that 
equations ( 2 ) yield: 


A = D = l/(l - V 2 ) 

B = v / (l - v 2 ) 

F = l/[2(l + V jTj 

y 

With sufficiently long members and a plate length L 
assumed that 


(3) 

it may be 

(4) 


w = Y 6 in(pitx/L) 


2Q 
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in which Y is a function of y only and p is the number of half 
waves in the X- direct ion. Insertion of equation (4) in equation (l) 
yields 

D - 2(B + 2F)X 2 + (AX 2 - 0 2 )\ 2 Y = 0 (5) 

dy 4 dy 2 

in which 


X - pjt/L 
0 2 = ta cr /EI 




(6) 


With Y = e a ^ this leads to 


w = (C ± cosh a^y + C 2 sinh a^y + cos a^y + Cj^ sin a^y) cos 


piTX 

a 


( 7 ) 


with 



( 8 ) 


in which 


q = B + 2F 
D 

H = (B + 2F) 2 - AD 



N 




(9) 


On the other hand for the restraining plates, with thickness t', 
a similar equation such as equation (l) also applies. Considering, for 
example, the lower half of a flange, as sketched in figure 5, this leads, 
analogous to equation (7)> to the general equation of the deflection 
surface 


xO 
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w' = (C 1 ' cosh o^'y’ + C 2 ' sinh o^'y' + C 3 ' cos Og'y’ + 

C^' sin ag'y ' )cos (10) 

with 

o^g' = j^GX 2 + X|/hX 2 + K(0') 2 (11) 

in which ( 0') 2 = t'a cr /EI' an< ^ 1 ' = (t')Vl2* 

Since the buckling and restraining plates are rigidly connected, 

X = pxt /L has, of course, to be the same for all plates. 

Dealing first with the restraining plates, for the sections dealt 
with in this report it is only necessary to consider the case of figure 5* 
Let the buckling plate exert a moment My’ = M cos(pnx/L) on the 
restraining plate. In the plastic domain the bending and twisting 
moments in the restraining plate are (references 3 and 4) 


V 


= -El' 


6 2 w* 6^w* 

A 22L. + B ° 


6x" 


(dyO : 


My' ^ -El 

' = -2EI'F 


B — — tc- + D 




6x 2 (dy') 2 

.P d 2 w’ 

dx dy' 


= "V 


( 12 ) 


Hence at y' = 0, where w' = 0 and d^'/dx 2 = 0, the boundary condi- 
tions are, from equations (12), 


and 




' = -EI 'D 


w' = 0 

(13a) 

6^' pnx 

(13b) 

(vF L 
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At y' = t>' (fig. 5) the plate is free, so that My' =0 and the so- 
called equivalent load Py f =0. Choosing the positive directions in 

the same way as in reference 6, figure l6o, page 295 * and using the 
equation for the transverse shear from reference 6, page 297 > 


p y' 


, _ = . 




+ . T V - 

3x by' 6x 


6My X ' 


6x 


or from equations (12) 


V = EI ’ 


d^w' i T _\ b\' 


(dy' ) 


m3 


+ (b + 4f) 


bxrby' 


m 


Hence from equations (12) and ( 14) at y' = h' the boundary condi- 
tions My' =0 and p ' =0 are 


b^ + d^2 = o 
6x 2 (6y') 2 


(13c) 


» -^4 ♦ (B + *F) - = 0 (13d) 

(6y>) 6x^y' 

Inserting equation (10) in the four boundary conditions (l3a) to (13d) 
yields the constants cp, Cp', cy, an< 3 ^4' i n equation (10), so 
that the angular rotation at y' = 0 can be computed from equation (10), 

|p - ( c 2' a l' + <V°2’)co» £2 

* (vV' 4 )’ (15) 

The value l/i = ^My 1 ) yi.Q^dw'/dy’ )y'=o may Pe called the s P r i n S 
constant of the restraining plate. Hence 

e = EIDty 

= E^)(6v , /^y , )/M > r , 


(16) 


which value will be used in the boundary conditions of the buckling plate, may be computed, 
yielding (reference 3> equation (56), p. 60) 


(£) : 


[M^)* - 

(a 2 ') 2 (q') 2 ] 

sinh sin a^'b' - cx-^ f cx^ T | 

[(q 1 ) 2 + (r 1 )^] 

cosh a-^'b’ cos a^'b 

* - 2a 1 , a2 , q , r' 

i 

( a l’) 2 + ( a 2* 

>31 

ja^* (r* ) 2 cosh a^'b* sin a^ 1 

b' - a^’b' cos ,_ b 



where 


q’ = (o ^) 2 - (Bk 2 /D) 

r' = (a-2 1 ) 2 + (BX2/ D ) [ (l8) 

X = pn/L 

If a.2 ' becomes imaginary, with a-^ 1 = a^" and 02* = ia-g"; 

/t\3 R 01 !") 2 ^") 2 + (02") 2 ^") 2 ] 8lnh a l" b ' sinh 02 %' - ^"ttg" [(q") g + (r") 2 ~[cosh c^'V cosh gg"b' - 2a 1 "a 2 "q"r" 

^ ' [(“l") 2 ' (“a") 2 ] [^l"( r ") 2 cosh c^'V sinh a 2 'V - a 2 ,, (q") 2 sinh c^'V cosh a 2 "b 


ro 


where 


q" = (a^') 2 “ BX 2 /d 
r" = -(a 2 ”) 2 + Bk 2 /D 


(l8a) 


If more than one restraining plate like that of figure 5 meets at y* - 0 the total spring 
constant 



(19) 
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so that for a symmetrical flange l/i^ = 2/\(r and \(r^ = VM or from 
equation (l6) 

e t = (1/2)0 (20) 

Considering now the buckling plate (fig. Mb)), the restraining 
plates are at the same time the "beams" by which the buckling plate is 
supported. At y = 0 (fig. Mb)) the difference between the transverse 
shearing forces acting on an element dx of the left beam will have to 
be in equilibrium with the equivalent load Py transmitted by the 
buckling plate and the resultant of the compressive forces A 1 a cr acting 
on the beam element dx. Here A-^ is the cross section of the left 
beam. The pertinent equilibrium condition is given by the equation 
(reference 6, p. 3^-6): 


■R d W _ _ A rr d^w 


( 21 ) 


in which is the flexural rigidity of the left beam about its major 

axis. At y = 0 


p y % 


^22 

dx 


_ ^ dMxy | 6My X 

dy dx dx 


or, from equations (12), 


P y = -El 


l d § + <b + *& 


( 22 ) 


so that from equation (21) one of the boundary conditions of the buckling 
plate at y = 0 is (reference 3 > equation (i), p. 6l) 


d^w 

Bl a? + EI 


D & +(B + to) a!fe 


d 2 w 

+ A-L0 cr ^2 - 0 


v 


The other boundary condition at y = 0 follows from equation (l5). 
Figures 4(b) and 5 show that for a one-sided flange = (My') y , _o 


Ik 
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and (dwA>Oy=o = (dw'/dy ' ) y « =0 > so that equation ( 15) becomes 


I - -Vih 


(23) 


From figure h(b) M , is a negative moment. By equation (l6) ^ may 

he expressed in terms of 0^> For a symmetrical flange 0 y is given 
by 0 t in equation (20), where 0 follows from equations (17) or (17a) . 

The subscript 1 refers to the left flange (fig. Ma)). Hence from 
equations (12) and (l6) equation (23) yields as second boundary condition 
for y = 0 (reference 3, equation (il), p. 6l) 


el- — * — \ - 0 

37 ' 1 \ D 3x 2 3y 2 / 


In the same way the boundary conditions for y — b become 


52 5? 


EI 


D ^ + (B + 4F)-% 
ox^oy 


Sy3 


5^w 

+ ^ 3 ? 


= 0 


<3w . (b 'd 2 v _ _ 


The subscript 2 refers to the right flange. Insertion of equation (7) 
in the four boundary conditions for the buckling plate leads to four 
linear homogeneous equations, which yield only values of the constants 
different from zero if the denominator determinant vanishes. This leads 
to the general buckling condition (reference 3, equation (6l), p. 6 2) 
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jo.-^ 2 js^S2 “ a 2 2r2 " (®l s 2 + ®2 s l) r2 + ®l®2 r ^ J " 

a 2 2 

|^l s 2 ” (®l s 2 + ®2 S 1 ) < 1 2 + + ®l®2 s l s 2 T ^|' s: ^ nB a l^ s ^- n ~ 

a l a 2 |^ s l s 2 “ (®i + 9 2 ) ( S 1 + 8 2 ) t ^ + 2 ( 0 1 S 2 + ®2 s l) < 3 r + 

20 102 < 3 2r2 ] cosl:1 a l b cos a 2 , ° + a 2 T j a l 2 S 1 + s 2 “ ( 0 1 + 0 2) r2 j “ 

^0^ + 9 2) s l s 2 + ®1®2( S 1 + s 2^ s ^ n ^ 1 a l^ cos a 2^ > + a i T '| a 2 2 || 

■\ 

(s x + 0 2) q2 ] + ( 6 1 + ^2) s l s 2 ' 0 l 0 2( s l + s 2) : 
s l s 2 + (®l s 2 + ®2 s l) qr + ®l®2^ 2r ^j = 0 


1 + s 2 - 


cosh a^b sin a^b + 


in which 


q = 


- (b\ 2 /d) 


r = a 2 2 + (b\ 2 /d) 


t = q + r 


= a 2 + a- 2 


1 2 


> 


= 2X./h\ 2 + K0 2 
%2 = f^( B l,2* 2 “ A l,2 ff cr) 


( 25 ) 


i6 
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For the elastic range A , B, D, and F are given by equations (3), so 
that from equations (9) 


G = 1 
H = 0 

K = 1 - v 2 


(26) 


The same buckling condition (24) obtains here, in which, consequently. 


a lj 2 — + x/td cr /N 


q = a-^ 2 _ v \2 


2 ,2 

r = + vX 


r = 2x/t a cr /N - 

s l,2 = V( Bl ^ 2> " 2 ' A l, 2 a cr) 


N = 


El 


1 - v 


while in equation (17) 


a l,2' = l / ±A - 2 + X it ' a cr /r 

q’ = (a ,') 2 - V X 2 

r ’ = (a 2 ') 2 + vX 2 


N' = 


El' 


1 - v c 


)> 


(27) 


(28) 


For an I- or H-section where both flanges are alike 6-^ = 6r> = 9 
and 82= 82=3 by which equation (24) transforms to 
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3Q 


(a-^u 2 _ ag 2 v 2 


a l^ a 2^ T ^ + 0 2s2T2 )sinh a^b sin agb - 


2a l a 2( w ^ - 20 st 2 ) cosh a-jb cos agb + 

2agT ^a-^2 u _ Qsv) sinh a^b cos agb + 

2ajj (ag 2 v + 0suj cosh a^b sin agb + 2a^agn^ = 0 

in which 

u = s - 0r 2 
v = s - 0 q 2 r 

n = s + 0qr 


(29) 


(30) 


while Q is given by 0-^ from equations (20) and (17). 

If the deflections of the beams (the flanges) in their own plane 
are neglected, so that s = oo, equation (29) reduces to (reference 3, 
pp. 59 and 63): 


a^ tanh ^o-jb/2) + ag tan^Xgb/2^ + 0t = 0 (31) 

This condition applies for any symmetrical rotational restraint of the 
web. In reference 3 the pertinent expressions for 0 were derived for 
several cases. For I- or H- sections 0 is given by 0+ from equa- 
tions (20) and (17). 

With relatively narrow and thin flanges, the rotational restraint 
exerted by the flanges on the web may be neglected, so that 0 = 00 and 
equation (29) transforms to 


( a l 2r ^ - Og 2 ^ ” s 2 -!- 2 ) sinh ajb sin agb - 2a^a 2 q 2 r 2 ^cosh a^b cos agb - l) + 

p p 

2a 2 q st sinh a-jb cos agb - 2a^r st cosh a-jb sin agb = 0 (32) 


For a T-section, where the web may be considered as the buckling 
plate, at the free edge of the web sg = 0 and 0g = m. Denoting s-^ 


l8 
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and 0p, which refer to the restrained edge, by s and 0, equation (24) 
transforms to (reference 2, equation ( 3 * 0 ) 

^a-j_ 2 r 2 u - a 2 2 q 2 v)sinh a^b sin a^o - a 1 a 2 (q 2 u + r 2 v)cosh a-|b cos a^b + 
a 2 T ^ a l 2r2 " ® sq2 ) sinil a i b cos + 

a l T ^ a 2 2 q 2 + 0sr 2 j cosh a^b sin a 2 b - 2a-ja,2qrn = 0 (33) 

where 0 is given by equations (20) and (17)* 

For an angle with equal legs the legs may be assumed to transfer no 
bending moments My to each other, so that 0 = °°. Consequently from 
equation (33) the buckling condition for each flange becomes (reference 7, 
equation (4o)): 


^a 1 2 r 1 *’ - a 2 2 q^ sinh a^b sin a 2 b - 2a 1 a 2 q 2 r 2 (cosh a^b cos a 2 b - l) + 

a 2 q 2 sT sinh a^b cos a-gb - a-j_r 2 sT cosh a^b sin a 2 b = 0 (34) 

In case of a T- stiffener or a column of a section as given by fig- 
ure 6, at edge 1 sj_ = °° and 0-^ = a>, so that, denoting S2 and 9 2 

at edge 2 as s and 0, equation (24) transforms to 
0 st 2 + a-^a2T 2 coth a^b cot a. 2 b + 

a 2 T (0 q 2 - s) cot a 2 b - a^T^er 2 - s) coth a^b = 0 (35) 


This case was completely worked out in reference 8, chapter 2.2. It 
was found that in case of a T-shaped steel sheet stiffener, which buckles 
at the yield stress, the interaction between the two modes of buckling, 
that is, torsional and local buckling, was negligible. 
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DERIVATION OF INTERACTION FORMULAS BY METHOD OF 
SPLIT RIGIDITIES 


The method of split rigidities divides the elastic or elasto -plastic 
behavior of a composite structure into its component parts. After the 
individual buckling stresses for these component parts are calculated, 
an interaction formula, which combines the individual critical stresses, 
gives the actual critical stress of the composite structure. 

This method was previously applied to the calculation of the critical 
stresses of built-up columns (references 8 and 9) and to flexural and 
torsional buckling of angles (references 7 , 8, and 10) and open sections 
in general (reference 10) as well as to the calculation of critical 
stresses in sandwich plates (references 1 and 11). The application to 
sandwich plates was very extensively explained in reference 1. In its 
application to the interaction of column and local buckling the method 
has to be used in a more generalized form. If two or more modes of 
deformation are involved in the buckling process, such as column buckling 
on the one hand and plate buckling of the web on the other in the case of 
an I-section, these two types, denoted as cases (1) and (2), respectively, 
are first considered separately. While one type is considered, the 
rigidity against the deformation of the other type is assumed temporarily 
to be infinite. 

In considering cases (l) or (2) separately, an equation may be 
established between the internal reactions and external actions. For 
case (l) above, for example, it is appropriate to compare the internal 
and external bending moments. In other cases, for example, case (2) 
above, it is convenient to establish an equation between the restraining 
and deflecting forces acting on a small element of the plate. In general 
for each case the most appropriate internal and external actions should 
be compared, which may, for other cases, differ from those mentioned 
above. For each separate case the external actions (bending moment, 
deflecting force, or otherwise) are directly proportional to the buckling 
stress for that case (exp or o 2 for cases (l) or (2), respectively) and 
to the deflection with buckling (w x or w g , respectively) . Hence these 
external actions may be expressed in these values: a lf Wj_ or a 2 , w 2 . 

Since during incipient buckling the internal actions are equal to the 
external ones, the former are also expressed in w^ and cr 2 , w 2 

for cases (l) and (2), respectively. 

For the actual combined case, with buckling stress cr cr and deflec- 
tion w = w^ + w 2 , the deflections w^ and w 2 are of form similar to 

those for the separate cases or are assumed to be so. Consequently the 
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internal actions caused by these deflections are assumed to be also 
equal to those for the separate cases. Hence the total internal action 
may be expressed in terms of crp, a 2 , w^, and w 2 . 

On the other hand the total external action with actual combined 
buckling is proportional to the actual buckling stress n cr and may be 
expressed in terms of cr cr and of a linear function of the deflections 
w 1 and w 2 . Writing down the equalities of internal and external 
actions, firstly for the particular actions considered in case (l) (in 
the chosen example the bending moments) and secondly for those considered 
in case (2) (here the restraining and deflecting forces acting on a 
small element) two homogeneous linear equations in w^ and W2 are 

obtained, from which the buckling condition is found by equating the 
denominator determinant of these equations to zero. Since this deter- 
minant contains o^, cr 2 , and o cr this gives a formula which expresses 
a cr in terms of the known individual buckling stresses cr-^ and a 2 * 

In order to show more directly how this method is applied, the 
explicit case of an I-section column will be considered first. In the 
section entitled "General Case of Columns with One or Two Planes of 
Symmetry" the general case of sections which are symmetrical with one or 
two axes will then be dealt with. In the sections entitled "Buckling of 
Tubes with Square Cross Section" and "Buckling of Columns with H-Shaped 
Cross Sections" cases of special sections, that is, tubes and H-sections 
which show some special features, are examined. It follows that for all 
these sections the interaction is in general negligible. 

The interaction is important only if the web that is perpendicular 
to the direction of column buckling is simply supported or elastically 
restrained at one side and free at the other side, so that a combination 
of flexural and torsional buckling occurs. This case occurs with 
T-section columns and with angles. These sections are dealt with in the 
sections entitled "Buckling of Columns with T-Sections" and "Buckling of 
Angles with Equal Legs." 


Buckling of I-Section Column in Direction Perpendicular 

to Plane of Web 

Derivation of interaction formula .- First an I-section with com- 
paratively narrow and thin flanges is considered (fig. 7(a)) so that, 
when the web buckles, practically no rotational restraints are exerted 
on it by the flanges. 

If the col umn buckles in a direction perpend icular to the plane of 
the web, its deformation may be split into two parts: (l) Buckling as a 
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column, without distortion of the cross section (fig. 7(t>)) and (2) dis- 
tortion of the cross section (fig. 7( c ))* 

If only case (l) occurs with buckling in a half wave of length a 
the buckling stress is 


a 


1 " 


JT 2E t 

( a /r) ^ 


( 36 ) 


where is the tangent modulus and r is the radius of gyration of 

the cross section. In the elastic range, = E. 

With a deformation according to case (2) alone, where, as stated 
above, the web is practically simply supported at the unloaded edges, 
for the same half wave length a the buckling stress would be (refer- 
ence equation (37) or reference 1, equation (3^-)) 


°2 = 


jt 2 EI 


b 2 t 


(a/3 2 ) + 2(B +2F) + DP^J 


( 37 ) 


or in the elastic domain 


Jt 2 N/l ' 
° 2 = + 3 , 


where 


) 

P = a/b 


( 38 ) 


( 39 ) 


In order to obtain the interaction formula, an equation will first 
be derived which expresses the equality between the external moments and 
the internal moments in the cross sections of the column. 

If only a deformation according to case (l) occurs, involving at an 
arbitrary point a deflection w^_ (figs. 7(a) and 7(t>)) it follows from 

the equality of internal and external moments denoted by and Me, 

respectively, that = M e or 


M il = p l w l = A c a l w l 


(40) 
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in which A c is the cross-sectional area of the column. This equation 
holds in all cases where the deformation w^ has the shape of a sine 
wave, no matter whether it is produced by an actual critical load Pj_ 

or by some other agent which causes the member to deflect in this shape. 
This is so because the internal moment is merely a function of the 
curvature. For the case of sine-wave deflection this curvature is 
proportional to the deflection itself at all points, so that internal 
moments are given by equation (4o), no matter what the cause of the 
sine -wave deflection. 

This will also hold true for the plastic range, if P]_ or a-]_ are 
calculated from equation (36) by using the tangent modulus corresponding 
to the actual axial stress a cr in the column. 

Assume now that a complete deformation according to figure j(a) 
occurs by buckling with a half wave length a. Then by the deformation 
according to case (l) the internal moment is given by equation (4o). By 
the deformation according to case (2) (fig. 7 ( a )) "the deflection of the 
web will vary in the Y-direction according to a curve which may be 
approximated sufficiently accurately by a sine wave. (This would actually 
hold true for w-j_ = 0.) Hence the average deflection of the web from 

case (2) is about (2 /ji)w 2. With concentric cylindrical buckling of 
the single web with a half wave length a the buckling stress for the 
elastic range would be 





( 4 l) 


or from equation (38) 



1 



( 42 ) 


NACA TN 261+0 


23 


Hence , in the same way as shown above , it follows from the equality of 
internal and external moments that, with a cylindrical deflection ^ 2 c> 
Mi = M e or 


M i = A v ct 2c w 2c 

= a w t 1 <j 2 w 2c (43) 

where A w is the cross-sectional area of the web. Therefore the 
average deflection (2/ir)w 2 from case (2) causes an internal moment 

M i2 = ! V ff 2 w 2 {kk) 


Consequently with buckling according to figure 7(a) the total internal 
moment is , from equations (1+0) and (1+1+) M-^ = + M^ 2 or 


Mi = A c 


a l v l 


2 

+ n r ^ 


From equation (1+2) 


n = 


l/(l + P 2 ) 2 


( 45 ) 


(46) 


With the additional notation 

2 K 

n A c 

Equation (45) becomes 


(47) 


Mi = M n +. M i2 


= A c (ct 1 w 1 + +I0a 2 w 2 ^ 


(48) 


If the actual critical stress is ° cr > deflection w^ from case (l) 
causes an external moment M e ^_ = A c cr cr w*L while deflection W 2 from 
case (2) causes an external moment M e 2 = A w a cr — W 2 * Hence the total 
external moment is 
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Me = M el + M e 2 

= k c °cr M l + Vcr § w 2 ^ 

or, using equation (47) , 

M e = A c o cr (w 1 + 0w 2 ) (50) 

Since Mg and Mi have to be equal, equations (48) and (50) yield, 
since A c cancels. 


a cr( w l + 0 W 2) = a l w l + T 10 a 2 w 2 (51) 

which gives one equation for finding the interaction formula. Another 
equation can be obtained by expressing the equality of the deflecting 
and restraining forces acting on a small element of the web. 

If only a deformation according to case (2) occurs (fig. 7(c)), 
that is, if w-l in figure 7(a) were zero, the critical stress in the 
elastic domain is a 2 from equation (38). For the equilibrium shape 

of the middle plane in case (2) at an arbitrary point P of the plate 
the deflecting force acting on a small element t dx dy is 
D 2 = -ter 2 dx dy(d 2 w/dx 2 ) . 

Since for the same shape of the middle plane at that particular 
point P the second derivative d 2 w/dx 2 is proportional to the maximum 
deflection w 2 of the plate, D 2 is proportional to a 2 and to v 2 
and hence it may be denoted by ca 2 w 2 . Here c is a proportionality 

factor which is a constant only for the given point P and for the same 
shape of the deflection surface, but is different for different points of 
the plate. Since the restraining force R 2 ac ting on the element and 
caused by the transverse shear stresses is equal to the deflecting 
force D2 it follows from the equation 

R 2 = 1 2 

that 

R 2 = cct 2 W2 (52) 

In the elastic domain the restraining force depends only on the 
shape of the middle plane of the plate. In the plastic range the ratio 
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of internal moment to curvature depends on the modulus and thus on the 
magnitude of the average compressive stress. Hence equation (52) remains 
true in the plastic domain even if the compressive stress differs from 
provided that in calculating a 2 from equation (37 ) > values A, B, 
D, and F are calculated for the actual compressive stress a cr . 

Considering again the actual deformation of the column during 
buckling according to figure 7(a) , it follows that the deflection wp 
from case (2) causes a restraining force Rp according to equation (52). 
The deflection w^ from case (l) causes a cylindrical bending of the 

web. If only this cylindrical bending occurs, the buckling stress of 
the web is equal to c 2c from equation (42). In that case the average 
deflection force acting on various elements of the web would be 

greater than ccr p c w^ because the proportionality factor c for a given 

point P of the web, as stated above, applies to a shape of the deflec- 
tion surface similar to that of w 2 (fig. 7(a)). The influence of the 
deflecting forces is the greater the greater the distance of the elements 
on which they act is from the edges of the web. Hence their influence 
may be expressed fairly well by comparing it with that of continuously 
distributed loads (proportional to w^ or w 2 ) on the bending moment 
in the middle of a simply supported beam with span b. For a uniform 
load q representing a distribution like w 1 , in the middle (y = 0 in 
fig. 7(a)) this bending moment is i qb 2 . It is A- qb^ for sinusoidal 

8 t[2 

distribution like Wg. Therefore as an average the deflecting force D-^ 
on an element in case only deflection w^ occurs may be expressed as 

bp = (n2/8)ca 2c Wj. With the notation 


7 = * 2 /8 (53) 

it follows then from the equality of restraining and deflecting forces, 
or R^ = D-^, that 

R ! = c 7cr 2c w;L (54) 

Consequently the total restraining force acting on an element of 
the web with buckling according to figure 7(a) is, from equations (52) 
and ( 54) , R = R 2 + R 2 or 


R = ca 2 w 2 + c 7 a 2 c w 1 


(55) 
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or, using equation (42), 


R = C0 2 (w 2 + T 17 w l) (56) 

With the actual critical stress a cr deflection w 2 causes a 
deflecting force co cr w 2 while deflection w 1 causes a deflecting 
force c 70 cr wi. Hence the total deflecting force is 

D = ccr cr (w 2 + Tw-^ (57) 

Equating this deflecting force to the restraining force from 
equation ( 56 ) yields, since c cancels, 

a cr( w 2 + 7w X ) = cr 2 (w 2 + t^) (58) 


Equations (51) and ( 58 ) are two equations with three unknowns, 
a CT , w x , and w 2 . This is necessary and sufficient, since w-^ and w 2 

have a common arbitrary factor, so that only the ratio w 2 /wj_ has to 
he known. Writing equations (51) and ( 58 ) as follows: 


(01 - 0 cr )w 1 - 0(0 cr - T)CT2) W 2 = 0 
7(0 cr - ri0 2 )w 1 - (a 2 - tf cr )w 2 = 0 


(59) 


the "buckling condition is obtained by equating the denominator deter- 
minant of these equations to zero, yielding 


J cr 2(1 - 70) ‘ 


0 n + (1 - 270 t i)0 2 - 


'J0 X + (1 - 2.7</>r\)a^ 2 - 4(1 - 70)(o x - 70R 2 a 2 ) a 2 
in which, from equations (46), (47), and (53) 


( 60 ) 


q = 


(1 


1 



2 


(46) 
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* - 1 % (k7 > 

n 2 

7 = g- (53) 


so that 


70 = -p — = 0.78 — 
^A c ' A c 


(61) 


while cr^ and are given by equations ( 36 ) and ( 37 ) or ( 38 ). 

Comparison with exact theory .- In order to check the accuracy of 
equation ( 60 ) its results will be compared with those of the exact 
calculation, according to which the buckling condition for this case 
is given by equation (32). It is given also, though in different nota- 
tions, by equation (t), page 3 V 7 of reference 6. For several channel 
sections with narrow flanges figure 182 , page 3^j of the same reference 
gives values of \|f against the ratios a/b = (3. With the notations of 
the present report 



A check of this graph shows, however, that the moment of inertia of the 
supporting beams, the flanges, with a width d and a thickness t, 
has been erroneously assumed for these channels to be equal to -i- td3 

12 

instead of to half of that of the entire cross section with respect to 
the axis of inertia parallel to the web. Hence the curves actually 
apply to I-section columns with flanges of width d rather than to 
channels. The curves for d = k inches and d = 2 inches are reproduced 
here in figure 8 , where also the pertinent cross section is given. The 
Euler curves for the same I-sections, plotting a ^ according to equa- 
tion ( 36 ) with E^ = E, are given in figure 8 by the dashed curves. All 
curves refer to the elastic range. 

Although the web width b should theoretically be measured between 
the middle planes of the flanges by which it is supported, for a correct 
comparison it will be measured here between the inner faces of the flange 
as was done in the exact calculation in reference 6 . According to equa- 
tions ( 36 ) and ( 38 ), with E+ = E and for a flange width d = 4", a 

Off u 

web width b = 15 ^ , a plate thickness t = 5/1 6", and Poisson's 
ratio v = l/k. 
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cJi/E = 4.50/a 2 

t 2 /e = 0.000345^ + 


while from equation (6l) 


70 = 0-52 


(63) 


(64) 


(65) 


For several ratios 3 = a/b values of q, o- j_/E, and cr 2 have been 
calculated from equations (46), (63); and (64), respectively. Hence 
a cr /E is computed from the approximate interaction equation (60) . The 

results for d = 4 inches are given in the table below. Corresponding 
approximate values \|r from equation (62) have been computed and inserted 
in the table as ty cr . These values have to be compared with the exact 
values -Jr according to figure 182 of reference 6 or figure 8 of this 
report, which values are also given in the table. 


d 

(in.) 

3 = — 

b 

T) 

W 

1 — 1 

b 

Q 

ro 

fcT 

a cr/E 

\|r cr 

^ex 

Percent 

error 

4 

1 

2 

2.5 
3 

3.5 
6 

0.25 

.04 

.019 

.01 

(.0057) 
( .0007) 

0.0181 
.00452 
.00289 
.00201 
.00148 
. 00050 

0.00138 

.00215 

.00290 

.00383 

.00495 

.01312 

0.001338 

.00166 

.001706 

.001522 

.001246 

.000469 

6.20 

6.90 

6.99 

6.60 

5.98 

3.67 

6.20 

6.90 

7.00 

6.70 

6.05 

3.72 

0 

0 

-.14 

-1.50 

-1.15 

-1.35 


It follows that equation (60), derived by the method of split 
rigidities, is very accurate and deviates but slightly to the safe side. 

For an I-section with d = 2 inches, according to equations (36) 
and (38), 

CT l/E = °-T3/a 2 (66) 

cr 2 / E = 0.000345^3 + jr) 2 (67) 

while from equation (6l) 


70 = 0.62 


( 68 ) 


MCA TN 2640 


29 


For this cross section the pertinent values have "been calculated in the 
same way and are given in the table below. Values of \|r ex have been 

obtained from figure 182 of reference 6 or figure 8 of this report. 


d 

(in.) 

3 = a 
b 



or 2 /E 

a cr/F 

^cr 

^ex 

Percent 

error 

2 

1 

0.25 

0.00294 

0.00138 

0.001147 

5.73 

5.85 

-2.2 


2 

.04 

.000735 

.00215 

.000611 

4.18 

4.08 

2.5 


4 

(.0035) 

.oooi84 

.00624 

.000176 

2.24 

2.24 

0 


Here, too, the agreement between \|r cr and \|r ex is very satisfactory. 

In figure 8 values \|/ cr from equation (6o) are given by circles showing 

the excellent agreement with the curves. The cross section of figure 8 
is unusual, since the moment of inertia with respect to the minor axis 
is very low as compared with the cross-sectional area. 


Simplified, more approximate interaction formula .- For normal 
I-sections or channels the values of a^, and cr cr , if plotted 

against the ratio 3 = a/b of half wave length a to web width b, 
will vary as in figure 1. The buckling stress is governed by the 
horizontal line AB for ratios 3 between Pj_ an< ^ $2 anc * ^y the 

curve for a cr for about 3 > Pg- At A, for 3 = 3]_, a-]_ is very 
high as compared with erg* For 3 = Pqj on "the other hand, is 

very high with regard to Equation (6o) may be written as 


a-, + (1 - 270 h) Cg 

/ 

1 

^ 4(1 - 70) (op - 70 ti 2 ct 2 )ct 2 

acr 2(1 - 70) 

-L ~ 

w 

/ p-L + (1 - 270 h)ct2 2 

y 


(69) 


If now ct-^ is much higher than such as is the case for 3 = Pp the 

last term under the radical is very much smaller than unity. Denoting 
this term by e> the term in the braces can be written as follows: 


1 






4* • • 


(70) 
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in which 


4(1 - 70)1 

- 70ti 2 ct 2 

Q 

ro 

1 

[a 1 + (1 - 270 ti)cx 2 ] 

2 


(71) 


From the derivation of equation (69) it follows that conservative values 
for a cr are obtained if t| is equated to zero, so that from 

equation ( 71 ) 


4(1 - 70)aicr 2 

(“1 ♦ * 2 ) 2 


(72) 


Hence from equations ( 69 )^ (70)^ and (72) 


o 


cr 


x °JL_ + U - 

°1 + “2 1 (<>1 + Or) d 


(73) 


Since » cr 2 , in the last term in the brackets a )^ j ( °1 + °2)^ ma y 
be equated to unity, so that equation ( 73 ) reduces to 


°cr ** ( 1 


70 


^1 + a 2 


o ± + (l - 70) cr 2 

a i + a 2 


(74) 


In the same way at P ^ (3^^ where as a rule 02 » it is 
found that 


cr 


1 - 70 


a l + a 2 / 1 


Or, + (1 - 70) cr-, 

= -*= 4 cr, 

rr_ 4- rr_ -L 


O 2 + cr 1 


( 75 ) 
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Since yj) is always smaller than unity, this yields for 3 = 3^ as 
well as for 3 > £9 the conservative values 


a l a 2. 

o = — — — — — 

cr a;L + a 2 


■ ("i' 1 + “2' 1 )' 1 (76) 

Criterion for range of negligible interaction . - From equation (74) 
for P = P 1 the decrease of cr cr with respect to cr 2 due to interac- 
tion is 


(An) 2 


a 2 " °cr 


7^2 
0 1 + a 2 


a 2 


(77) 


From equation (75) for P = P 2 the decrease of o cr with respect to 
crp is 


(Act) 2 = a 1 - a CT 
= 1 — a, 

a 2 + a i 1 


(78) 


At 3=32 value cr 2 is equal to Oj for 3 = P 2 (fig* l), so that 
from equation (77) 


(Acr) g 


70(0i)p = p 2 

( a l)p=P 1 + (-l) p= p 2 02 


(79) 


From equation (36), since P = a/b 


(<Tl) P=Pl 



(80) 
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It should be noted that this equation, as well as equation (36), 
neglects the influence of shear deformations, as is customary. However, 
for such small buckling lengths as considered here (for 3 = 3p the 
buckling length a »b) the shear deformations may substantially decrease 
the buckling stress cr^ • It was shown in an unpublished paper that this 

does not influence the resulting governing interaction as expressed in 
equation (93)* From equations (79) and (80) 


(A a) 2 = 


70 


(P2/Pl)‘ 


+ 1 


(81) 


On the other hand, at 3 = 3 2 the value is equal to cr 2 for 

3 = 3p (fig. 1), so that from equation (78) 


(Aa ) 1 


^( g 2)|3 = p 1 

(° 2) P={ , 2 + (° 2 ) Wl 


From equation (38) 


(° 2 ) 


3=32 



) + 3 2 

Liv^) 

+ 3 X _ 


(02) 3=3i 


(82) 


(83) 


so that from equation (82) 


(Aa) x = 


7$_ 


"( 1 /P 2 ) + P2 

2 

+ 1 

_(i/pi) + fi 



(84) 


The stress cr 2 from equation (38) is minimum for 3 = 3p = 1- Assuming, 
for an example, the rather low ratio 32/Pp = 10, so that 3g = 10, and 
using equation (6l), equation (8l) for 3 = 3^ yields 


(Act) 2 


0.78A w /A c 
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Equation (84), for p = P 2 , gives 



= 0.030 J- cr 1 


( 86 ) 


L c 


With A w /A c =l/3 or 2/3, equations ( 85 ) and (86) give decreases 
(A 0)2 and (Acr)^ of 0. 0026(72 or 0*005202 and 0.010-^ or 0.020^, 
respectively, and thus maximums of l/2 and 2 percent, respectively. 

From equations (8l) and (84) the interaction is the smaller the greater 
the ratio 32/3p, so that for P2/Pl = 10 it is negligible. In the 
elastic range 32 is determined by the condition (fig. l): 


From equation ( 38 ), since is minimum for 3 = 3p = 1, with v = 0-3 


^ ai ^3=32 ( a2 )3=3x 


(87) 




( 88 ) 


so that from equations ( 87 ) and ( 36 ) 



or 


( a / r )p = p 2 = 1 . 65 b/t 


(90) 
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or 


P2/Pl = 



1.65r/t 


(91) 


Hence, since for p 2 /Pl - 10 the interaction is negligible, it is also 
negligible for 


Approximate formula for maximum amount of interaction .- If P 2 /Pl 
is of the order of magnitude of 10, it is seen from equations (8l) 
and ( 84 ) that for other ratios P2/P1 than 10 the discrepancies (Aa) 2 
and (Aa) 1 will vary practically proportionally to (P1/P2) 2 ° r > from 

equation (91), to (t/r) 2 . From equations ( 85 ), (86), and (92) for 
Pg/P}. = 10 or r/t = 6 the governing discrepancy Aa is 3(A V /Ac) p< 

cent. Hence in general the maxim-urn decrease of or a 2 hy intera- 

tion will be 


It was assumed above that the rotational restraints offered by the 
flanges to the web could be neglected. It may be shown, however, that 
equation (93) remains valid if the web is substantially rotationally . 
restrained. From the derivation of equations (51) and ( 58 ) it is evident 
that these equations and hence the resulting buckling condition (60) and 
its further elaboration apply for that case as well. Only 0, 7, and q 

change in value in that case. But even equation (6l) for 70 remains 
practically the same. Indeed, the restraints caiise a decrease of 0 
but an increase of 7, so that a calculation shows 70 to be practically 
independent of the rotational restraints. Even for negative restraints 
(see the section entitled "Buckling of Columns with H-Shaped Cross 
Section") a,nd also for the quite different case of the web of a T-section 
(see section entitled "Buckling of Columns with H-Shaped Cross Section ) 

70 from equations (ll8) and (17^) is O.QA^/P^ and 0 . r J^>A- w /A C} respec- 
tively, and thus practically equal to 70 = 0 . 78 A w /Ac from equation (6l) . 
In calculating the individual buckling stress a 2 of the webs the 



(92) 



(93) 
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restraints at their unloaded edges have, of course, to be taken into 
account. In this case q will have a value different from that of 
equation (46). However, q may always sufficiently accurately be 
equated to zero, as was done from equation ( 72 ) on. 

Finally equation (93) is based on equations ( 86 ) and (92). To see 
whether equation ( 92 ) remains valid for rotationally restrained web 
edges the most extreme case of built-in web will be considered. For 
this case (reference 6) the minimum value of Og is 


jt% 


(02) 3=3l = 7 $£ 


b2 t 


■ 

and occurs for a/b = = 0*7* Hence now, instead of equation ( 89 ) 


_(a/r)2 

or 

(a/r) p= p 2 = 1.25b/t 
or 

P 2 = ( a / b )p=p 2 = l-25r/t (95) 

With r/t =6 equation (95) yields P 2 = 7*5, so that 

^ 2 /Pl = 7‘5/0. 7 = 10 . 7j and hence it is but slightly higher than 

^ 2/^1 = 10 for r/t = 6 and simply supported web edges. 

Furthermore, for a given ratio P 2 /P 1 , (A cr) 2 /a 2 from equations (79) 

and ( 36 ) is independent of the rotational restraints of the web. On the 
other hand (Acr)i/aj from equation ( 82 ) will decrease with increasing 
restraint. From reference 12 for fully clamped edges 



0 2 = [^/3 2 ) + 2.5 + 5P^^ 


(96) 
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Hence for r/t = 6, P 2 = 7-5 and Pp = 0.7* and, using equation (6l), 
equation (82) yields 


(Aer)p 


0.78A w /A c 

41.6 


0.019 



(97) 


Consequently from equations (86) and (97) for r/t = 6 the governing 
discrepancy A a = (Aa)-^ for varying amounts of web restraint will vary 

between 3(A W /A C ) percent and 1.9(A W /A C ) percent, so that for any posi- 
tive rotation restraints equation (93) yields safe results. Apparently 
with negative web restraints, such as occur in H-sections, for r/t = 6, 
Act may be larger than 3(A W /A C ) percent, so that equation (93) errs 

on the low side. This will also become apparent from the section 
entitled "Buckling of Columns with H-Shaped Cross Section." 

If cr cr is in the plastic range, for simply supported web edges 
from equations (87), (36), and (37) one has, instead of equation (89), 


(a/r f 




or, instead of equation ( 91) , 


E t r 

P 2 = 1.65 zr- J 

qE X 


(98) 


Here Ep is the tangent modulus and q is the reduction coefficient 
for plastic plate buckling (references 3 and 4), both referring to the 
actual buckling stress cr cr . Since Ep/E is smaller than q, for 

equal ratios r/t the value of P 2 is smaller than in the elastic 
range. However, P p is likewise smaller than its value Pp = 1 in 
the elastic domain (reference 3). Nevertheless, P 2 /Pp may be somewhat 
smaller than in the elastic domain. 

Equations (79) and (82) apply also to the plastic domain, where crp 
and cr 2 are given by equations (36) and (37)> respectively. Equa- 
tions (80) and (8l) also remain valid so that for the same ratio P2/Pp 
the ratio (Acx) 2 /ct 2 for P = Pp is the same as in the elastic domain. 
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Because in equation (37) factors (B + 2F) and D are always higher 
than factor A, the ratio between the values of cr^ for (3 = P 2 

and (3 = 0^, for the same ratio P2/PIJ is even higher than that 
expressed by equation ( 83 ) for the elastic domain. Hence for equal 
ratio P 2 /P 1 the ratio (AcOp/cr^ from equation ( 82 ) is still smaller 

than in the elastic domain. 

Moreover all discrepancies between a cr and cr^ or °2 are 
greatly reduced in the plastic range because with decreasing cr cr the 

pertinent plastic values E-^, A, B, D, and F increase. This may 
reduce (Acr)^ from equation (93) many times, in the same way as a 
relative difference in elastic buckling stress may be reduced many times 
in the plastic range. For steel columns or plates of which the buckling 
stress is, for example, at the yield stress (in the sense that the 
tangent modulus is zero) an increase in the elastic buckling stress by 
100 percent will result in no increase in the actual plastic buckling 
stress. 

Therefore in the plastic domain equation (93) will yield safe 
value s for ( Aa ) ^ . 

General Case of Columns with One or Two Planes of Symmetry 

From the preceding derivations for columns with an I-section it is 
evident that it makes no difference whether the symmetrical support of 
the web is effected by symmetrical flanges or otherwise. Therefore 
formula ( 93 ) applies in general for any column of which the pertinent 
web is symmetrically supported and positively rotationally restrained 
at both unloaded edges and which buckles in the plane of symmetry and 
in a direction perpendicular to that web. Examples are columns with 
an I-section, channels, and box sections. In the latter case (fig* 9)> 
which will be dealt with more extensively in the section entitled 
T, Buckling of Tubes with Square Cross Section" in the case of buckling 
in the mode of figure 10, A w refers to the joint cross section of the 
two horizontal webs, which buckle individually in the same direction as 
the entire column. 

Moreover, from the derivation of the interaction formulas (7*0j 
( 75 )> and ( 7 6 ) it follows that they apply as well for columns which 
buckle in the plane of the web, and which are symmetrical with respect 
to that web. Examples of this kind are H- and T-sections. The plate 
buckling stress in the above-mentioned formulas refers here to the 

buckling stress of a plate which is clamped at one unloaded side and 
free at the other. For this case 70 follows from equations (157)> 
(l6o), (l 66 ), ( 167 ), and (175) and is equal to 0.69Af/A c where Af is 
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the total cross section of the flanges. For these cases too the interac- 
tion between column and plate buckling is generally negligible. 

The interaction is important only if the pertinent web, that is, 
that perpendicular to the direction of column buckling, is elastically 
restrained at one side and free at the other, so that a combination of 
flexural and torsional buckling occurs. This case is dealt with in the 
section entitled "Buckling of Columns with T-Section." 


In the square tube the four plates act as if they are simply supported 
at the unloaded sides, so that their buckling stress in the elastic domain, 
from equations ( 38 ) and (88), is 


In this case condition (92) is certainly satisfied so that the 
interaction will be very small. However, in connection with the tests 
on these sections reported later herein, they will be studied in some 
more detail. 

The following two cases will be discussed separately: (a) The 

plates buckle at about erg. One has then to investigate whether such 
plate buckling induces column deflections and consequent interaction, 
(b) The member buckles by column deflection at about crp, that is, at 
P 2? pg in figure 1. In this case one has to investigate whether such 
column buckling induces plate deflections and, if so, what their effect 
is on interaction. 

(a) The plates buckle symmetrically with respect to the vertical 
and horizontal axes of inertia of the tube (fig. 11), so that the total 
equivalent load p y which they transfer to the column as a whole is 

zero, in the vertical as well as in the horizontal direction. Hence 
no bending of the column as a whole is induced by local buckling of the 
plates. It is as if the column were infinitely rigid against bending, 
so that in figure 1 at p = P]_ = 1 the situation is the same as if 
were infinite. Hence at P = P^ in figure 1 a cr is not smaller than 
a 2 so that 


Buckling of Tubes with Square Cross Section 



( 88 ) 



( 99 ) 
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and no interaction occurs. (Of course, a very slight deflection of the 
unloaded edges of, for example, the horizontal plates occurs because 
the vertical plates in figure 11 will compress very slightly in the 
transverse (vertical) direction by the loads transferred to them by the 
buckling (horizontal) plates, but the resulting deformations may be 
neglected. ) 

At 3 = p 2 (fig- 1) likewise the buckling stress of the 

column is not diminished by a buckling of the cross section according to 
figure 11 because this does not change the position of the center of 
gravity of the cross section. Hence it does not cause an extra external 
moment. Neither does it cause a change of the internal moment. Con- 
sequently it does not affect the buckling stress of the column, so that 
at 3 = 3 2 


ff cr “ CT 1 “ 


ir 2 E 


(a/r)‘ 


( 100 ) 


and no interaction occurs. 

Interaction may occur, though, with buckling of the webs according 
to figure 12. Here, however, cr 2 is much higher than according to 
equation (88) . The plate buckling stress in this case is the same as 
for a section according to figure 12 with b 1 = b/2. The buckling 
stress for this case was calculated (reference 3) from equation (3l)> 
where, as was derived in reference 3 


6 = e a 



coth a^b* - a 2 f cot a 2 f b f 
( a l') 2 + ( a 2') 2 


and a-^ T and a^ ! are given by equation (ll) above. The ratio b/t 

at which the horizontal plate buckles at a given stress is expressed in 
terms of those ratios (b/t)^ and (b/t) ss plates with both 

sides fixed or both sides simply supported, respectively, would buckle 
at the same stress. This relation is expressed in reference 3 by 
equation (62) 



( 101 ) 


tb * 

where 7-, for several ratios p. = -—rr- is given in figure 13 of that 
-L t b 

reference. It follows from equation (88) that 
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3.62 ^ = 1.90 



(102) 


For a plate with both unloaded sides fixed from equation (9^*) 


For the section of figure 12, with t* = t and b f = b/2, 

[i = tb f /t f b = 0.5, whence from figure 13 of reference 3 one finds 
7l = 0.52. Thus equation (105) yields 


Since the factor 7 ^ for the horizontal plate in figure 12 is about 

midway between those for simply supported (7p = 1) and for fully 
clamped plates ( 7 ^ = 0), the optimum half wave length will be about 
the same as for a plate with one side clamped and the other simply 
supported, namely, 0.80b. Hence here P]_ in figure 1 is 0.80. 

For a square tube (fig. 9) the moment of inertia and the cross- 



(103) 


Hence from equation (101) 



(104) 


or 



(105) 



(106) 


sectional area are 



P P 

and A c = B - A , respectively 


that the radius of inertia is 



r 


12 


(107) 
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To take an example, with B = 2.5 inches and t = 0.062 inch, 

A = 2.376 inches and r 2 = 0.991 square inch. With 
a = 0.80b = 0.8 o(B - t) = 0.80 X 2.438 inches = 1.952 inches, equa- 
tion (36) yields 


CT 1 = 


(l.952) 2 /o.99l 


E 


2.56E 


While from equation (106) 


a 2 = 4.84^|^|) 2 E = 0.00314E 

Since equations (7^-) and (76) from the section entitled "Buckling of 
I-Section in Direction Perpendicular to Plane of Web" apply in the present 
case, from the more conservative equation (76) at P> = Pj_ 

°r-r = 0.00314E = 0.003136E (108) 

cr 2.563 


It is apparent that the interaction is negligible even for the mode of 
figure 12. Moreover, since for the mode of figure 11 cr cr from equa- 
tion (99) is smaller than for that of figure . 12, buckling according to 
figure 12 will not occur in the first place and a cr is given by 
equation (99), resulting in 


a cr - 



0.00235E 


(109) 


(b) With buckling as a column at 0 = 02 the deflection wj of 

the column causes deflective forces -tcrc) 2 w/8x2 (fig- 10) which cause 
a similar deformation of the cross section as in figure 12. For an 
equally distributed vertical load q on the horizontal plates of fig- 
ure 10 the moments M in the corners C follov from the equation 


qb^ Mb [ Mb f _ 2 Mb 
2 51 ~ 2N 3N ~ 3 N 


yielding 
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or 3/A times the clamping moments for a fully clamped plate. Hence a 
conservative value for the buckling stress is obtained by using the 
average of a simply supported and a fully clamped plate. Thus from 
equations (38) and (96) 

a 2 = |l/(B 2 ) + 2.25 + 33^^ (HO) 

At p = (3 2 the value of is equal to ffg for 3 = Pi so from 

equations (88), (89), and (90) 



and 


3 = a/b = 1.65 | (111) 

As A = b - t and B = b + t, from equation (107), since t/b is small, 
r = O.Ulb, so that from equation (ill) 

p 2 = O.678 | (112) 


or with the dimensions given above 


- °- 678 ■ 26 - 6 


Hence from equation (110), for 3 = 32 


cx 2 = 2129.25 


2.A38 2 x 0.062 


N = 1.2414-E 


In figure 1 at 3 = f rom equations (88) and (109) 

*1 = ( a 2) p=p = 0-00235E 


that from equation (7 6) a conservative value for a cr at 3 = 32 


°cr 


1 . 2kb 

1.24635 


0.00235E = 0.002316E 


SO 


is 
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Hence the effect of interaction is here less than 0.2 percent. According 
to the more accurate equation (75)> since 70 is here about 
0.75A W /A C = 0.75 x 0.5 = 0.38 

a cr = 0.00235E = 0.002348E 

Hence the more accurate effect of the interaction is only 0.1 percent. 

The approximate formula (93) yields an effect of 

l (o. 4 ib/t ) 2 percent = 2 (1/0^062 ) percent = °- 2 percent 

This formula gives values which are too high because the mode which 
yields the minimum value c^, that is, plate buckling according to 
figure 11, does not govern the interaction in this case of buckling 
according to figure 10. 

Thus for concentric buckling the effect of the interaction between 
column and local buckling in square tubes may be neglected. 


Buckling of Columns with H-Shaped Cross Section 

Somewhat more -pronounced interaction between column and plate 
buckling may occur in H-sections where the width of the flanges is such 
that they are rotationally restrained from buckling by the web (fig. 13)* 
In that case the optimum half wave length of buckling will be larger than 
the width b of the web, so that in figure 1 p^ > 1. Thus the 
ratio P 2 /P 1 on the interaction depends may be smaller here than 

if the web were simply supported or elastically rotationally restrained 
by the flanges and therefore more interaction may occur. 

Let the H-section have the dimensions given in figure 13, so that 
b = A - t = 2.875 inches, b ! = 0.5B = 2.3125 inches, and t = l/8 inch. 
The buckling condition is given by equation (29), but it may more easily 
be obtained by using equations ( 7 ^-) and (75), which, according to the 
section entitled "Buckling of I-Section Column in Direction Perpendicular 
to Plane of Web," apply here. 

The plate buckling stress a 2 follows from reference 13 from the 
condition 


L S = S IV + 2S 111 = 0 
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which, using the pertinent tables (reference 14), is satisfied for a 

minimum stress 


a 2 = 0.002008E = 21,500 psi 


(113) 


for a half wave length a = 2b = 5-75 inches. Thus here = 2 (fig. l) . 
The radius of gyration of this section with respect to the middle plane 
of the web (fig. 13 ) is given by 


2 x 

r 


ItB 3 


-L At 3 
12 


1 4.4 


(2B + A - 2t)t 


1 

6 2B 


B' 


+ A - 2t 


2 ( c T^) \ ^ 

or numerically r^ = 1.375 square inches, so that ^ = » * 

l* 3 1 

Hence from equation ( 36 ), with E^. = E, at P = (fig- l) 


2k. 


a-, = - * .--5 = 0.411E (ll4) 

1 (a/r) 2 

The same value 70 as found in equation (6l) for the case of figure 7 
will approximately apply here , so that 


70 * 0.78 


A - t 

2B + A - 2t 


0.185 


Hence equation (74) yields 


cr 


cr 


0.411 + 0.815 X 0.002 

07413 


0 2 


= 0 Q ^^ 3 0.002008E 


so that at 


= 0.002006E 

(3 = Pi the interaction is negligible. 


( 115 ) 


= 0.002008E 


At p = (3 2 


CT 1 = 


n 2 E 


2 ( CT 2)p = p. 


(a/r)^ ' “'^1 

so that, with r 2 = 1.375 square inches, from equation (ll6) 

= /4920r 2 = 82.3 in. 


( 116 ) 


Hence 


a 


p 2 -r - = 28.6 

2 h 2.875 


, _ a _ 82.3 _ 


p * - ^ - 

2 h' 2.3125 


= 35.6 


Since for these high (3 values tables are not available, the plate buckling stress cr 2 
for (3 = p g ) and hence for the case where the deflections of the beams are neglected, was 
calculated directly from the pertinent equation (31) 

tahh^a^b/2^ + a g tan^a 2 b/2) + 0r = 0 (31) 

in which 0 is given by 0^. from equations (20) and ( 17 ), 


1 / t ^ ( a x' 2r ' 2 ~ a 2 ,2q ' 2 ) slnh a i' b ' 3in a 2 'h' - a-p'o-g' (^' 2 + r ' 2 )cosh o^'b' cos a^'b 1 - 2a 1 l a 2 l q , r t 
2 * ( a l' 2 + a 2 ,2 )( a i' r ' 2cos h a-p'b' sin a^'V - a 2 'q ,2 sinh cos ag'b'j 

f 

By trial and error it was found that equation (31) is satisfied for Og = 0.0675E. 


VJl 
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Moreover, this case may also be calculated by the method of split 
rigidities. This method yields a buckling stress of O.0665E. Assuming 
thus 


a 2 = O.0675E (117) 

the actual critical stress can be calculated from equation (75)* Since 
with the assumed large half wave length in the X-direction the web is 
bent approximately according to a parabola by the moments My exerted 

on it by the more unstable flanges, it follows from the derivation of 
equations (47) and (53) that here 

0 =f A w /A c = 0.16 


7 = 1.2 


70 = 0.8 A w /A c = 0.192 
Hence from equations (75 ) > ( ll6) , and ( 117 ) > 


cr 2 + (1 - 70)a-, 

cr cr = cr-, 

+ a n x 


0.0675 + (0.808) (0.002) „ 
0.0695 1 


(118) 


= 0 . 001996 E 


Even here the interaction diminishes only by 0.5 percent. Since 

here Aw/A c =0.24 and r/t = 01. 375/0. 125 = 9*4 ^a rough estimate from 

equation (93) would yield a decrease of 3 — (— ^ percent = 0.3 percent. 

Ac \r/t / 

That equation in the section entitled "Buckling of I-Section in Direction 
Perpendicular to Plane of Web" was derived for webs with no restraints 
or with positive rotational restraints only. It was already remarked 
there that for negative web restraints it will give values for (A o) ■$_ 
which are too small. In the present case, however, the order of magnitude 
is correctly indicated by equation (93)* 
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Buckling of Columns with T-Section 

In the stability of such sections it is best to distinguish two 
cases: (a) The flange has a relatively small torsional rigidity or is 

not significantly more stable than the web (fig. 14) or, more generally, 
the web is not substantially rotationally restrained by the flange, and 
(b) the flange has a relatively large torsional rigidity and is much 
more stable than the web, so that the web is substantially rotationally 
restrained. 

On the basis of a numerical example it will first be shown that in 
case (a) the interaction from cross-sectional distortion is practically 
nil. Case (b) will then be discussed in greater detail. 

(a) If the ratio b/t of the web is not much greater than the 
ratio b f /t ! of the flange and the flange is relatively weak (fig. 14), 
so that the web is only slightly rotationally restrained by the flange, 
both web and flange remain practically straight in cross section; that 
is, the cross section is not distorted during buckling. Hence in that 
case a column with T-section buckles approximately in the way that it 
is assumed to buckle in the analysis of flexural and torsional buckling. 
The interaction formula for the buckling stress for this case was derived 
by Kappus (reference 15) from three simultaneous different equations of 
the fourth order and also by Lundquist and Fligg (reference l6) . It was 
also derived by the first author according to this method of split 
rigidities (reference 10) . The exact solution for this problem was 
given by the first author in references 2 and 3> including the influence 
of distortion of section, while in the present report it is given by 
equation (33)« 

In reference 10 the buckling stress for a steel column with an 
effective length L of 700 centimeters and with the cross section of 
figure 15 vas calculated from the interaction formula of flexural and 
torsional buckling (without regard to cross-sectional distortion) as 
well as from the exact equation (33). The corresponding buckling stresses 
were 9^8 and 942 kilograms per square centimeter, respectively. (in 
reference 10 the first stress is also given as 9^2 kg/cm 2 , which value 
was due to a computational error.) Hence in this case the exact buckling 
stress is only 2/3 percent less than that found from the assumption that 
the cross sections of the plates do not bend. 

It was stated in reference 10 that this is due to the high value of 
the critical stress of a plate which is clamped at one unloaded side and 
free at the other for a half wave length of many times its width, as 
compared with the buckling stress of the T-section. The additional 
influence of this plate buckling stress can be taken into account in the 
following manner: The torsional -flexural buckling stress is firs ' t 
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computed according to the interaction formula for this case. If, in 
analogy to the derivation of equation (7 6) in the section entitled 
"Buckling of I-Section in Direction Perpendicular to Plane of Web," 
this mode with o 1 _ 2 is regarded as case (l), case (2) is represented 
by the above buckling of the web clamped at one unloaded side and free 
at the other, with a corresponding buckling stress This is true 

because such buckling is the only deformation which could occur if the 
column were infinitely rigid against the deformation of case (l). The 
transverse bending of the flange may be neglected. Hence, as will be 
shown in detail later, an interaction formula results which is again 
approximated conservatively by a formula similar to equation ( 76 ). The 
actual combined buckling stress is then given by 


0 cr 


q 3 

cr 3 + 


a l-2 


(119) 


in which a j_-2 = 948 kilograms per square centimeter = 13,460 psi and 
is the buckling stress of a plate that is clamped at one unloaded 
side and free at the other (reference 12). 


°3 " 



+ 0.125P 2 


+ 



(120) 


where (3 = a/b = L/b. In this example the width b of the web, which is 
actually supported in the middle plane of the flange, is 30.5 centimeters, 
(fig. 15), so that (3 = 700/30.5 = 22.9, yielding, with v = 0.3, 


a 3 " (527 + 0,57 + 




b^t 


= 66.472 


12(l - v 2 )b 2 t 


= 0.0645E 


= 0.0645 x 2,100,000 kg/cm 2 = 135,500 kg/cm 2 


= 1,925,000 psi 
Hence from equation (119) 

°cr - 948 Wc* 2 = 941 k«/cm2 = 13,360 psi 
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which checks the exact buckling stress of 94-2 kilograms per square 
centimeter = 13*375 psi found from equation (33). 

It is seen that for this particular example: (l) The interaction 

formula, equation (119)* is accurate within about l/lO percent, and 
(2) the interaction between torsional-flexural and plate buckling 
amounts to about 2/3 percent, that is, is negligible. 

(b) If the flange has a relatively large torsional rigidity and is 
much more stable than the web, as, for example, for the dimensions of 
figure 16 , the web is rather strongly constrained rotationally by the 
flange. In this case the influence of plate buckling, with a buckling 
stress Oy is more pronounced, especially for small half wave lengths. 

At the extreme, with a flange which is infinitely rigid, only pure plate 
buckling could occur. 


The deflection of the cross section consists here of a translation wp, 
a rotation with respect to the shear center S, denoted by the deflec- 
tion w 2 at the lower side of the web and a plate buckling of the web 
with a maximum deflection w^ (fig. 17) • The shear center may be 
assumed sufficiently accurately to be situated at the intersection of 
the middle planes of flange and web. Again, the transverse bending of 
the flange is negligible. 


If only the deflection occurs, that is, flexural buckling only, 

the critical stress is 


or in the elastic domain 


rt 2 E. 


CT-, = 


( a / r y) ' 


a l = 


i t 2 e 


( a / r y)' 


( 121 ) 


( 122 ) 


This case will be called case (l). If only the deflection wg occurs 
(case (2)), the critical stress is (reference 10) 


Oo = J-(: 


EFI+ + E. C 


P\ 




or in the elastic domain (reference 17 ) 


(123) 


(124) 
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where I p , GI^., and EC W are the polar moment of inertia about S, 
the torsional rigidity, and the warping rigidity, respectively. Equa- 
tion (123) follows directly from equation (124) since in the plastic 
domain, if during buckling infinitesimal shear stresses T xy' are 

superimposed on pure compressive stresses, the ratio r xy' between 

excess shear stresses and excess shear strains, which is G in the elastic 
domain, is replaced by EF (reference 3). Furthermore, in equation (124) 
the elastic modulus E refers to excess bending stresses, so that in the 
plastic domain it changes to the tangent modulus E t . 

The deflection of the web with a maximum w^ (case (3)) would, if the 
column were infinitely rigid with respect to the other deflections, 
correspond to a critical stress in the plastic domain: 


In the elastic range equation (125) reduces to equation (120). The 
critical stress in the plastic domain according to equation (125) may 
be found by the energy method of equating the internal work Vj_ to the 
external work V e , hence 



(126) 


in which (reference 3) 



2 


+ 2B 



and 



(128) 


Assuming 


v = w^IjL - cos (jty/2bj] sin (jtx/a) 
equations (126), (127), and (128) yield 


(129) 



(130) 
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This gives for the elastic domain, where, from equations (3), 

A = D = 1/(1 - v2), B = V/(1 - v2) # and F = l/[2(l + v )] , 
with v =0.3, 

a 3 = ^§”^5 + + 0 . 68 ^ ( 131 ) 

instead of the more -accur ate value from equation (120). This is due to 
the fact that the energy method yields values which are too high. Hence 
equation (130) may be improved by adjusting it to the mofe -accurate 
known value in the elastic range. This is done by multiplying the coef- 
ficients of D by 0.125/0.138 and those of B and F by 0. 57/0.68. 

This is how equation (125) was obtained. 

Some years ago the first author derived the interaction formula for 
the critical stress a cr as expressed in terms of the critical 
stresses cfp, ag, and a ^ for the component modes. This was done by 

considering successively (l) the equilibrium between external and internal 
bending moments, (2) that between external and internal torsional moments, 
for the entire column, and (3) the equilibrium between external and 
internal torsional moments acting on the web alone. This leads to a 
system of three linear homogeneous equations in Wp, Wg, and w^, from 

which the buckling condition follows by equating the determinant of the 
system to zero, yielding a cubic equation in o cr * 

In order to check this equation, during the summer of 1948 Mr. C. D. 
Maussart, the first author’s assistant at the Institute of Technology at 
Delft, Holland, calculated cr cr exactly for the elastic range from equa- 
tion (33)> where dp and 0.2 are given by equations (27), for the 
T-section of figure 16 with b = 31 centimeters, b* = 15 centimeters, 
t = 1 centimeter, and t* = 2 centimeters. The freely supported length 
was assumed to be 700 centimeters and a cr was calculated for elastic 

buckling in 1, 2, 4, 8, and 12 half waves and hence for half wave 
lengths a of 700, 350, 175> 87. 5> and 58.33 centimeters, respectively. 
The exact critical stresses a cr in pounds per square inch thus found 
are plotted in figure 3 against the half wave length a. In table 1 they 
are given in kilograms per square centimeter as well as in pounds per 
square inch under the column heading "a cr from equation (33)." 

The stresses <Jp, cr 2 , and a 3 from equations (122), (124), and 
(120) for the component cases (l), (2 ), and (3) are likewise given in 
figure 3 and table 1. The above-mentioned cubic interaction formula 
in cr cr was found to yield results which deviate less than 1 percent 
from the exact values. This proved the accuracy of the method, but did 
not yet provide a simple and direct way of computation, since the coef- 
ficients of the cubic equation are rather involved. 
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In the following discussion simpler interaction formulas are derived 
for this type of buckling and are checked against the accurate values from 
equation (33)- 

In order to find explicit interaction formulas for the specific case 
of p = (3-j^ (fig- l), the critical stress °2-3 first be calculated 

for v-l = 0, that is, for fixed shear-center axis, so that a 2 _ 3 follows 

directly from a quadratic equation like equation (60) in the section 
entitled "Buckling of I-Section in Direction Perpendicular to Plane of 
Web." Since cr^ here very large (in fig. 3 for a = 60 cm), so that 

a will be only slightly smaller than a 2 _o, it will subsequently be 
cr 

sufficiently accurate to express the interaction of a 2-3 an< ^ a i i n 

a simple formula of the type of equation (j 6 ) . On the other hand, at 
P = p 2 (fig. 1), hence, in figure 3 at a * 600 centimeters, w 3 is 

first assumed to be zero, so that the critical stress ^ 2.-2 ^ or fl exura l 
and torsional buckling is obtained as interaction of cases (l) and (2). 
Subsequently the small interaction between °i _ 2 and the very high plate 

buckling stress cr^ ma y again be expressed by a formula of the type of 
equation (76). 

Critical stress for T-section with fixed shear-center axis . - Although 
actually the shear center is situated slightly below the middle plane of 
the fl ang e it is sufficiently accurate to assume the shear-center axis at 
the intersection of the middle planes of flange and web. 

In a similar way to that done in the section entitled "Buckling of 
I -Section in Direction Perpendicular to Plane of Web" the entire section 
is considered first. Since, however, twisting about the shear axis occurs 
here instead of free bending, no equilibrium between bending moments can 
be considered. The best procedure is to compare the moments about the 
shear-center axis S of the deflecting and restraining forces acting on 
a small slice of length dx of the column (fig. 18). 

Assume first a deformation according to case (2) to occur which 
corresponds to w 2 - At buckling the deflecting force -tcr 2 (8 2 w/dx 2 )dx dy 
acting on a small element t dx dy of the slice in the web is propor- 
tional to o 2 and to w, so that it may be denoted by ca 2 wt dy. Since 
t dy represents a small element dA c of the cross section of the column, 
the deflecting forces on elements in the web as well as the flange may be 
denoted by cOgW dA c . 

Hence, since from figure l8 for case (2) the displacement w of any 
point of the cross section is £ w^, if r i s if s distance from S, 
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the total moment about S of the deflecting forces acting on the slice 
A c dx is 


M 


d2 


ca^wr dA 


c 


£ 

h 



dA 


c 


= — a w 0 I 
h 2 2 p 


(132) 


where I is the polar moment of inertia about S. 

P 

Consequently the moment M r about S of the restraining forces 

acting on the slice, that is, of the transverse shearing stresses acting 
in its end cross sections, is likewise 


M. 


r2 


= - CTpWpI 

b d d P 


(133) 


Assume now a deformation according to case (3) alone (corresponding 
to w^, fig. 18 ) or, in other words, assume the column to be infinitely 

rigid against a deformation of case (2). Since here the web alone 
deforms, if this deformation were of the same shape as in case (2) the 
deflecting moment about S would be given in equation (132), where Ip 
would be replaced by the polar moment of inertia Ip W of the web 

about S and c^w 2 by a 3 v 3* However, since with equal deflection 

at y = b the deflection at other points is less in case (3) than in 
case (2), the deflecting moment is also less, so that here 


M 


d3 


= a b ®3 W 3V 


( 134 ) 


where a < 1. The restraining forces are exerted here by shearing 
stresses in both end sections of the slice but also by the bending 
moment My exerted by the flange on the web at its built-in unloaded 


edge. This is an internal moment so that it does not add to the 
restraining forces acting on the slice as a whole. In equation (120) 
for cr , My is accounted for by the factor 0.125P^* Omitting it for 
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the reasons just stated, the restraining moment for this case is 


•rj'"; 13 V3“3 


in which, from equation (120), 


( 

!l/p 2 J 

1 + 0.57 

(1/P 2 ! 

) + 0. 

125P 2 +0.57 


(135) 


(136) 


Hence, with both cases (2) and (3) occurring simultaneously, the 
restraining moment about S is 

Mr = Mr2 + **r 3 


' a 2 V 2 + ar] 3 0 3 V 3J b h 


(137) 


At the moment of buckling, with a buckling stress cr 2 _ 3 , 
deflecting moment about S follows from equations (132) and (134) by- 
replacing a 2 and cr^ by °2-3> 90 


M. = 


M d2 + M d 3 


/ P w \c T 

= a 2- 3 ( w 2 + a T7 V 3yb b 


(138) 


Hence from equations (137) and (138) the condition of buckling M d = My 
is given by the equation 


J 2-3 w 2 


Wo + a 


w 3 = a 2 w 2 + a — ^°3 V 3 
P / P 


(139) 


In order to obtain a second equation in w^ and w^ the equilibrium 

of a slice bt dx of the web alone is considered by comparing again the 
moments of deflecting and restraining forces about S. 
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If a deformation according to case (2) alone 
stress of the web is, analogous to equation (124) 


occurs, the critical 
for the entire section, 



+ EC 


ww 



( 1 * 0 ) 


where the second subscript w refers to the web. The deflecting 
moment, using the same proportionality factor as in equation ( 132 ), 
would be 


M d2v = F a 2w w 2 I pw (l4l) 

which follows directly from equation ( 132 ) by replacing a 2 and I p 
by cf 2 W and Ip W * Hence the restraining moment in this case is likewise 


M. 


L r2w ^ a 2w w 2^pw 


( 1 ^ 2 ) 


With a deformation according to case (3) alone the deflecting moment 
is already given by equation (13*0 • Since at this stage the equilibrium 
of the web is considered, the restraining moment My exerted by the 
flange must also be taken into account, so that now' the moment about S 
of the restraining forces acting on a slice bt dx of the web is equal 
to i n equation ( 13 *+) 


M r3w = 


a — cr 0 w 0 I 
b 3 3 pw 


(1*3) 


whence the total restraining moment 

^rw = ^r2w + ^3^ 

= ( a 2w w 2 + aa 3 w 3 )& V (l44) 

At the moment of buckling, with a buckling stress cr 2 o, the deflecting 
moment about S follows from equations (142) and (l43) by replacing 
a 2w and by <?2-3> so that 
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M dw " M d2 + M d3 


CT 2-3( 


w 0 + aw 0 ) — I. 


3J b P w 


(145) 


Hence the condition of buckling 


M. = M 
dw rw 


is here 


a 2-3( w 2 + aw 3) = a 2w w 2 + aa 3 w 3 (146) 

Collecting terms in w^ and w^, equations (139) an h (l46) become 

(&2 - a 2 - 3) v 2 + i ^ 3°3 “ cr 2 - 3)^ aw 3 = 0 
( a 2w - a 2-3) w 2 + ( a 3 - *2-3) aw 3 = 0 


(147) 


where 


^ - ^pw/lp 


(148) 


The buckling condition is obtained by equating the determinant of 
equations (lVf) to zero, from which 


2-3 2(1 - Ilf) 


a 2 + a 3 ‘ * ( a 2w + ^3 a 3) " 


^p 2 + a 3 - t(°2w + 2 - 4(1 - i|r)ff 3 (a 2 - 


(149) 


where for the elastic range a^, a °2 m } ^3* and ^ are Si ven by 

equations (124), (120), (l4o), (136), and (l4£), respectively. The 
stress a 2-3 '* ias ^ een calculated for the T-section of figure 16 for 

several half wave lengths a. The result is given in table 1 and 
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figure 3* It follows clearly from figure 3 that , approaches o ^ 

for larger wave lengths, where cr^ is very high. 


Buckling stress of column with T-section for p » (fig* 1).- 

In order to find the real buckling stress cr cr the influence of a 
deflection w-^ of the column (fig. 17) has still to be considered. 

Since for P = 0^ (fig. l) cr^ is very high, this influence is rather 

small. The stress cr cr will be calculated from the interaction of 
case (l) with deflection w.^ (fig. 17) and the combined case (2-3), 
considered in the previous section, with a buckling strefes a ^ ^ and 

a deflection Wg_^ (f^S* 19). 


This case is similar to that of the I-beam in the section entitled 
"Buckling of I-Section in Direction Perpendicular to Plane of Web." 
Considering first the equilibrium between external and internal moments 
analogous to equation (40) the internal moment from the deflection w^ 

from case (l) is here 


M il = A c°l w l 


(150) 


Analogous to equations (43) and (44) the internal moment from the 
deflection v 2-3 of the web would be found here as 


M i(2-3) _ Vl CT 2-3^ w 2-3 


(151) 


in which 


T) = (152) 

2-3 

where a is given by equation (4l). The average deflection of the 
web from case (2-3) is pw 2 _ y However, this internal moment ^(2-3) 

is somewhat ambiguous. In case (2), where a pure rotation of the section 
with respect to the shear center occurs, according to the definition of 
the shear center no internal bending moments will originate. But equa- 
tion (151) indicates internal moments for this case. This is a consequence 
of the fact that the center of shear S was assumed at the intersection 
of the middle planes of flange and web, while in reality it is slightly 
lower. Therefore it is better and safe to neglect here, the 

more so since it is extremely small in comparison with from equa- 

tion (150). Hence 
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(153) 


With actual combined buckling the external moment resulting from w-|_ is 


From w 2-3 analogous to equation (151) 


To find |i it is noted that for a deflection of the web from case (2), 
where the cross section remains straight, the average deflection 
Hw 2 = 0 . 5wg, so that p 2 = 0-5- For case (3) the cross-sectional distor- 
tion of the web may sufficiently accurately be approximated by that of 
a cantilever beam with uniform load q (fig. 20). Thus 


M e(2-3) " A w CT cr^ w 2_3 


(155) 


M = - - q(b - y) 2 


2 



M dy dy 



so that w 3 = (w) y=b = t 4 and 


w = y 4 - 4by 3 + 6b 2 y_ 2 w 

3 ^ 


3 


(156) 


Hence 



(157) 


bw 


3 


Thus in equation (155) 


0.5 > h > 0.4 


( 158 ) 
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From equations (154) and (155) the total external moment is 


NL = M„-i + M. 


el + w e( 2 - 3 ) 


= A c a cr f w i + *2-3] 


or with the notation 


0 = M-A w /A c 


= 


= A c a cr( w l + ^ w 2 - 3 ) 
Hence from equations (153) and (l 6 l), since Mg = 


(159) 

(160) 

(161) 


a cr( w l + 0 w 2-3) = CT i v i 


(162) 


Next, in the same way as in the section entitled "Buckling of 
I-Section in Direction Perpendicular to Plane of Web," a second relation 
is obtained by considering deflecting and restraining forces acting on 
an element t dx dy = dA w dx of the web. If only case (2-3) occurs the 
deflecting force may be written as 80 that, analogous to 

equation ( 52 ), the restraining force in that case is also 


R 2 = ccr 2-3 w 2-3 (163) 

If only deflection w^_ occurs, the deflecting force is c 7 a 2c w l' 

where 7 accounts for the fact that with equal maximum deflection at 
y = b, a larger average deflecting force results with w = w^ = Constant 
than with a variation of w with y as in case (2-3) • A good measure 
of this influence is the moment of the deflecting forces with respect 
to S. As a matter of fact, this means that the moments about S of the 
deflecting and restraining forces acting on a slice bt dx of the web 
are considered and equated. Denoting the deflecting forces on an 
element dA w dx by kcrw dA w , the deflecting moment is for case (l) 

M^ = kcr r Wjy dAy. 

= 0.54b^tow-, 


(164) 
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For case (2), where w = ^ W2, 


M 2 = ka 


wy dA 


w 


= 0.333kb 2 tcrw 2 (165) 

For case ( 3 ), where w is approximately given by equation (156), 


rb 


= ka 


vy dAy. 


= 0.289kb^taw 3 

Hence from equations (l64), ( 165) , and (l66) 


(166) 


or 


-Oil. >7 

0.289 0.333 


1-73 >7 >1-5 


(167) 


Thus analogous to equation (54), the restraining force R x = D x or 

R-j^ = c7a 2c w 1 (168) 

so that the total restraining force is, from equations (163) and (l68), 

R = R 2 + Rj or 

R = ccr 2-3 w 2-3 + c7a 2c w l (!69) 


or from equation (152) 


R = ca 2 _ 3 (w 2 _ 3 + qrwj (170) 

With actual combined buckling the deflecting force is ca cr w 2-3 
for case (2-3) and c7a cr w 1 for case (l), so that its total is 


D — ca cr ^w 2 _ 3 + 


(171) 
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Hence from the equality of D and R 


a cr( w 2-3 + 7w l) = a 2 - 3 ( w 2-3 + ^l) 


(172) 


Collecting terms with w^ and w 2_3 i- n equations (162) and (172), they 
become 


(ai - a cr )w 1 - 0a cr w 2 _ 3 = 0 

7 (a cr - ^2-3) W 1 - (a 2 _ 3 - a cr )v 2 _ 3 = 0 




> 




(173) 


These equations are nearly identical with equations (59)- The only 
difference is that here w 2-3 occurs instead of and that in the 

first of equations (173) ^ has been equated to zero because ^(2-3) 
from equation (151) is negligible. From equations (173) 


°cr = 2(l 1 - 70){ CT 1 + (1 " 70Tl)0 2-3 " 

\J^ 1 + (1 - r 0 n)CT 2 - 3] 2 “ ^ " r0)cricr 2 _3 j- 


(174) 


Here t] is given by equations (152), (4l), and (l49), while from 
equations (158), ( l 6 o) , and (167) one can only conclude 


, A w 

°*75 — >70 >0.69 — 


(175) 


To decide which value 70 to use it i.s noted that cr cr is the smaller 
the larger 7 0. (The quantity 70 is an interaction factor, which is 
equal, for example, to unity in the case of sandwich plates) . Hence a 
safe and only slightly conservative value is 


70 = 0.75 p- (176) 

A c 

In the same way as was stated in elaborating equation (60), con- 
servative and sufficiently accurate values of a cr are obtained if T] 
is equated to zero. With r\ = 0 equation (17^) is identical in form 
with equation (60), if Oq is replaced by a 2-3 - Hence, for the 
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particular case of 3 = 3y, that is, for very large (fig. l) equa- 

tion ( 17 ^) may be written, in analogy to equation (7*0 > as 


<*1 + (1 - 70) a 2-3 
0 cr = a ± + a 2 _ 3 °2-3 


(177) 


or more simply and conservatively, analogous to equation (76), as 


cr 


“l + «2-3 2 ' 3 


(178) 


Here is given by equation (122), a 2 ^ by equation (l49) and 

70 by equation (176). For the particular section of figure l6 

70 = 0.75 = 0.25. (Comparing equations (177) and (178) it follows 

that the sandwich-plate formula (reference 1) is indeed obtained if the 
interaction factor 70 = 1.) 


In table 1, a cr is computed for the T-section of figure 16 from 
equation (177) as well as from equation (178) and given under the 
headings "o cr from equation (177)" and "° C r from equation (178)," 
respectively. Comparing these values of cr cr with the exact values 
obtained by using equation (33)j it is seen that the governing value 
of a cr , that is, for 3 = 3^ (fig. l), corresponding here to a value a 

of about 58.33 centimeters, may be computed more accurately from the 
simpler formula (178). The results from equation (178) for a = 58.33 
and a = 62 centimeters are indicated in figure 3 by circles, and are 
seen to fit almost exactly the exact curve for a cr . For a = 87*5 
and a = 135 centimeters, the more-accurate equation (177) has to be 
used, since here the difference between a cr and ^2-3 becomes more 
important. However, in these cases a cr for buckling in one half wave 
has no practical value, since buckling in more than one half wave results 
in smaller buckling stresses. The resulting values are indicated in 
figure 3 by crosses. The agreement with the exact curve is excellent. 


Computation of critical stress a 1 _ 2 for flexural and torsional 
buckling of T-section .- Although the more general derivation of a 1 _ 2 
according to the method of split rigidities was given in reference 10, 
the pertinent critical stress for the special case of a T-section will 
be derived here directly. In a computation of flexural and torsional 
buckling it is assumed that no distortion of the cross section takes 
place. Hence only deformations according to cases (l) and (2) occur, 
while in an approximate sense the column is supposed to be infinitely 
rigid against a deformation from case (3). 


NACA TN 26k0 


63 


The bending with respect to the Y-axis is considered first. If 
only case (l) occurred, a deflection w^ would cause an external 
moment A c 0 i w i> so that this deflection entails an internal moment 

M il = A c a l w l (^9) 

Since it is a property of the shear center S that a rotation with 
respect to S yields a torsional moment only and no flexural moments, 
the internal flexural moment from case (2) is zero. Hence the total 
internal flexural moment is 


M i = M n = A c a l w l (l80) 

The external flexural moment is found by multiplying the critical 
thrust ^ or ^l exura l an -d torsional buckling by the displace- 

ment CC f of the center of gravity C . From figure 21 this is equal 

y Q 

to w^ + — Wp, so that 


Me = A c a 1 _ 2 fw ± + ^ w 2 ) 
Equating Mg and Mj_ yields the condition 

0 l-2 W 1 + IT W 2 ) = °l w l 


(181) 


( 182 ) 


Torsion with respect to the center of shear S is next considered. 
Here a slice of length dx of the column is considered, as under the 
section entitled "Critical Stress for T-Section with Fixed Shear -Center 
Axis," where it was found in equation (133) that for case (2) the moment 
with respect to S of the restraining forces acting on the slice dx 
can be represented by the formula 

M r2 = F a 2 w 2 I p 

Again, from the definition of the center of shear it follows that the 
deflection w^ from case (l) causes no internal torsional moments, so 
that the resisting moment as above defined, for the combined cases (l) 
and (2), is 
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“r = ”r2 * § Wp 


(183) 


From equations (132) or (183) it follows directly that at the critical 
stress 2 the deflecting forces from the deflection w 2 cause a 

moment 


M d2 ~ b a l-2 w 2 I p 


(184) 


From the second member of equation (132) the deflecting moment from w^ 
is found as 


M dl = I ca l-2 wr ^c 


= ca 1 _ 2 w 1 A c y 0 

so that the total deflecting moment is 


(185) 


M d = M dl + M d2 


h 


- co l. 2 IVo w l * b *2 


Equating M d and yields 


r 2 
r o 


or 


a i-2 y o W l + — w 2 = *2 — W 2 


/by G \ 

a l-2^ W 1 + W 2j = CT 2 w 2 


Collecting terms in equations (182) and (187), 

( a l - a l-2) w l - V a l _ 2 w 2 = 0 


by 

-| a 1 _ 2 v 1 - (a 2 - a 1 _ 2 )w 2 = 0 


(186) 


(187) 


(188) 
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from which 


1-2 


i]T - (y 0 2/r 0 2)] [°1 + "2 ' A°1 ' "a) 2 + SV« 2 /’. ! _ 


(189) 


where y Q is the distance between the centers of gravity and shear and 
r© is the polar radius of gyration about S, while cr^ and o ^ are 

given by equations (122) and (124). For the T-section of figure l6 
°l-2 given in table 1 and by the curve a ± 2 in & ure 3* 

Buckling stress of column with T-section for P = Pp (fig- !) •- 

For 3 = the critical stress a 3 for case (3) is very high (fig. 3), 
so that the exact critical stress a cr will be only slightly lower than 
°1_2 • Hence it may be expected that the sandwich formula, where 
100-percent interaction is involved, provides a sufficiently accurate 
and safe approximation, so that ®cr ma y 4>e calculated from the formula 


°3 

° CT ~ a 3 + a x _ 2 a l-2 


(119) 


Indeed for a = 700 centimeters, where from equation (33) the exact value 
of the buckling stress is 1880 kilograms per square centimeter = 26,700 psi, 
equation (119) yields a cr = 1883 kilograms per square centimeter 
= 26,740 psi. .This value as well as that for a = 6 20 centimeters is 
aiso given in table 1, while in figure 3 both values are indicated by 
circles and fit the theoretical curve of cr cr exactly. 

To check the method of calculation in general, a cr will be 
calculated more accurately. Considering first the bending moments in 
the column with respect to the Y-axis, according to equation (153) in. the 
subsection entitled "Buckling Stress of Column with T-Section for 3 * Pj 
(fig. l)" the total internal moment is 


MjL = A c CT l w l (190) 

The total external moment by deformation according to cases (l) and (2) 
is, analogous to equation (l8l). 


M e(l-2) 


^c°cr 




(191) 
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In equation (157) i n the subsection mentioned above it was found that 
in case (3) the average deflection of the web is ^ 



0. 4w^ 


so that for case (3) 


M e3 = Vcr (0 * 4)w 3 

yielding a total external moment 


(192) 


M_ = 


M e(l-2) + M e3 


A a 
c cr 


' y Q . Av \ 

W 1 + — w 2 + °* 4 aT w 3 


(193) 


In order to limit to two the number of unknown deflections and thus the 
necessary number of equations and the degree of the buckling condition, 
it is observed that from equation (l82) of the preceding subsection 


y 0 °1 

w n + -t— w 0 = — w, 

1 b 2 ^i_2 1 


so that, with the notations 


a = 


g l 

7 l-2 


K 

6 = r 

n c 


equation (193) becomes 

M g = A c a cr (aw 1 + 0.46w^ 
Hence from equations (190) and (196) 

CT cr^ aw 1 + O.40w^ = a-^Wj 


(194) 


(195) 


(196) 


( 197 ) 


NACA TN 2640 


67 


Next a slice bt dx of the web is considered. From equations (l64), 
( 165 ), a^id (l 66 ) of the subsection entitled "Buckling stress of a col umn 
with T-section for (3 tv (fig. l)" the moments about S of the 
deflecting forces acting on the slice may be represented as 


= 0.5kb 2 tcrw 1 


(164) 


Mg = 0.3331*^ ov 2 


( 165 ) 


M 3 = 0.289kb 2 taw 3 ( 166 ) 

for cases (l), (2), and (3), respectively. The internal moment origi- 
nating from a deflection w^ of the web is negligible for the large 
half wave length at {3 ^ If only w 2 occurs, the buckling stress 

is 0 ^ from equation (l4o), while with w 3 alone, it is a 3 from equa- 
tion (120). These stresses result in deflecting moments, which, from 
equations ( 165 ) and (l66), may be denoted as 


M d2w = 0*33300^2 (198) 

and 

M d3w = 0.289ca 3 w 3 (199) 

respectively. Hence the total restraining moment is 


^W M r2w + M r3w 

^d2w + ^d3w 
or 


M rw = c (°- 3330^2 + O. 2890 3 w 3 ) (200) 

Similarly from equations (l64), ( 165 ), and (l66) the total deflecting 
moment with actual buckling may be represented as 


M 


dw 


C0 cr ( 0 . 5 ^ + 0 . 333^2 + 0 . 289 W ) 


( 201 ) 
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so that from equations (200) and (201) 

(0.5W-L + 0.333*2 + 0.289 * 3 )tf cr = 0.333^*2 + 0.28903*3 (202) 


From equations (194) and (195) 

w 2 = (a - 1 )— w-l (203) 

y o 

On the other hand, the position of the center of gravity C (fig. 21) 
is determined hy the equation 

Ay | = A cy 0 

so that from equation (195) 

3^. = 2 — = - (204) 

^0 ^ 6 

and from equation (203) 


^ a - 1 

W 2 - 2 —e~ V 1 


(205) 


Hence from equation (202) 


cr 


^0.5 + 0.667 a - — ) *]_ + 0*289*3 


= O.667 a 1 Ogw*! + O.28903W3 


(206) 


Collecting terms in and equations (197) an -d (206) become 


(°! ~ acJ cr) w l " °'^ a cr v 3 = ° 

jl.33(a - l)(J2w - [ e + !*33(a - l)] ff C r] w l + 0.5789 (a 3 - a cr ) V3 = 0 


> 


( 207 ) 
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from which 


°cr 2(1.33 + 0.11a - G) ( a l-2 + a 3) “ 1# 33( a - 1 ) a 2w 



jl . 44a | 

( a l-2 + *3) 

| - 1.33(a. - 1)0^] 

2 - 5.760-^(1.33 + 0.11a - e)J 


(208) 


in which a and G are given by equation (195). Using, in equation (70) 
the first term only, equation (208) simplifies to 


cr 


j 

a 3 + a l-2 ■ °'925cT2 W [i - (<*i_2/ ff i)J 


J l-2 


(209) 


which results conservatively in equation (119) as given before; that is, 

ct 3 


J cr 


2 - 2 


(119) 


The stress a cr has been calculated from equations (208) and (119) 
for a = 350, 620, and J00 centimeters, respectively, and is denoted in 
table 1 by "cr cr from equation (208)" and "<? cr from equation (119)," 
respectively. It is seen that for a = 620 and 700. centimeters, where 
0 cr is the governing stress, equation (119) yields even more exact 
results than equation (208). Its values are indicated in figure 3 by 
circles. For a = 350 centimeters, where the difference between cr 1 _ 2 
and a cr is much larger, equation (208) gives accurate results, while 
equation (119) is somewhat too conservative, as could be expected. The 
value for a cr from equation (208) for a = 350 centimeters is indicated 
by a cross in figure 3* 


Buckling of Angles with Equal Legs 

For reasons of symmetry the legs do not rotationally restrain each 
other, so that the exact buckling condition for this case is given by 
equation (34). On the other hand the buckling stress for flexural and 
torsional buckling (reference 10) is given by equation (189). To illustrate 
the use of the interaction formula, the exact and the approximate buckling 
stresses will be computed for a specific example. 

Assuming a steel angle with a length L = 40 inches, a width 
b = 2 inches, and a thickness t = 0.1 inch, b/t = 20, so that it 
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buckles in the elastic domain. The flexural rigidity B of a flange 
against bending in its own plane, for use in equation ( 27 ), is here 

Eb^t/3 instead of Eb^t/l2, since the axial strain e x from the deforma 

tion by bending, which results in the deflection w^ (fig. 22), is here 
zero for the fibers at the shear center axis S. By trial and error 
equation (3*0 yields the exact buckling stress cr cr = 0.000875E (refer- 
ences 7 a nd 8) . 

In the approximate equation (l89), as applied to this case, cr p is 
the critical stress for flexural buckling about the Y-axis (fig. 22), 
since with torsional buckling of the section the center of gravity 
displaces in the direction of the Z-axis, resulting in bending about 
the Y-axis. The stress cr 2 is given by equation (l2^). 

Since the flanges are not rotationally restrained, they remain 
practically straight in cross section, so that, in contradistinction to 
the case of the T-section in the section entitled "Buckling of Columns 
with T-Section,” case (3) does not occur here. 

For an angle, y Q = b/(2\/2) and I p = 2tb3y/3, so that 

r Q 2 = Ip^(2bt) = b ^/3 and y 0 2 / r 0 2 = 3/8- Hence equation ( 189 ) for 
the case of an angle reduces to 


cr 


= 0.8 


+ Or- 


~ \l(°l - ff 2) : 




( 210 ) 


With Iy = ( l/l2) ( t\(2) (b\(2) 3 = (l/3)b^t and r ^ = (l/6)b 2 , from equa- 
tion (122) for the particular angle quoted above, with a = L, 


From equation (124), with Lj. 


p = 0. 004llE 

( a / r y) 

= — bt^, while the last term is negligible, 


= G ^ = A(|j = °- 0009625E 


\2 


so that from equation (210) a cr = 0.000879E. This is in close agree- 
ment with the exact result cr cr = 0.000875E from equation (34). 
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Buckling according to the considered mode governs only if it yields 
smaller critical stresses than for buckling of the angle with respect to 
its minor axis. For the latter mode the interaction is negligible and 
the pertinent moment of inertia is 




12 


— b^t 
12 


and hence Iy/4, 
Since this value 
buckling, a cr = 


so that the column buckling stress is cr^/4 = 0.00103E. 


is higher than that previously computed for torsional 

0 . 000879E governs . 


THEORETICAL RESULTS 


Comparison between the exact and approximate methods of calculations 
leads to the conclusion that the interaction between column and local 
buckling may be computed very accurately by the method of split rigidities 
in the following manner: 


(a) For columns where the web plates are supported at both unloaded 
sides, such as columns with I-, H-, or box sections or with channels, for 
buckling in a direction perpendicular to the web or webs the interaction 
between column and local buckling is always practically negligible. As 
a way of estimating roughly the effect of the interaction, formula (93) 

may be used, which states that o QT is at most about 3 (a w /a c ) j^6/(r/t 2 

percent less than the column or plate buckling stress, o^_ or ( a 2) m i n > 

whichever is smaller. The total cross section of the web or webs located 
perpendicular to the direction of buckling is A^.. The total cross 

section of the column is denoted by A c . The radius of gyration referring 
to the above-mentioned direction of buckling is r, while t is the web 
thickness. 


The critical stress a Qr may be calculated from equations (7^-) 
and (75) more directly. For P = (3^ (fig. l) 


a cr 


+ (l - 70)cr 2 

a ^ + C2 


(74) 


where and o 2 are the column and plate buckling stresses, respec- 

tively, for the considered half wave length of buckling. This half wave 
length has to be chosen such as to make a cr a minimum. It can practi- 
cally always be taken equal to the half wave length that makes O 2 
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minimum. The value (a cr ) . thus obtained governs as long as p is 
smaller than about P 2 (fig* l) • 

For 3 > P 2 


0 

cr 


a 2 + (1 - 70)^ 

02 + CT-^ 


a l 


(75) 


where hoth a-^ and cr 2 have to he calculated for a half wave length 
equal to the effective length of the column. The value of 70 is about 
0-75A w /A c . More conservatively, in general 


'cr 


q l g 2 
CT 1 + °2 


= (a^ 1 + a^ 1 )" 1 ( 76 ) 

Equations (7*0> (75)j and (76) obtain also for the plastic domain, 
if the plastic values E-t/E, A, B, D, and F from equations ( 36 ) 
and (2) are calculated for the actual buckling stress a cr . Moreover, 
equations (7*0.? (75)> and ( 76 ) apply for buckling in the plane of the 
web for columns which are symmetrical with respect to the web, as for 
example H- and T-sections. In that case 0 2 refers to the plate buckling 

stress of the flanges which, for the considered direction of buckling, 
behave as plates which are free at one unloaded side and clamped at the 
other. In this case 70 is about 0-7Af/A c , where Af is the total 
cross section of the flanges. 

In all above-mentioned cases the interaction is practically always 
negligible . 

(b) The interaction between column and plate buckling is important, 
however, if the pertinent plate subtends a significant angle to the 
direction of buckling and is elastically rotationally restrained at one 
side and free at the other side, such as in T-sections and angles. In 
this case a combination of column buckling and twisting occurs. This 
follows directly from the conservative formula ( 76 ), because in this 
case both the flexural and torsional buckling stresses and 02 

have their smallest value for the largest possible half wave length, 
that is, the effective length of the column. For many practical sections 
these two stresses are of the same order of magnitude. For example, 
with 0 ^ = 02 , equation ( 76 ), though too conservative in this case, yields 
a cr = 0.5<?i = 0 . 5 ^ 2 ^ which shows how significant interaction can become. 
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If the plates do not restrain each other significantly rotationally, 
as in T-sections with rather weak flanges or in angles, the interaction 
is rather accurately accounted for by the formulas for flexural (column) 
and torsional (plate) buckling. Then 



where and are the column and torsional buckling stresses, 

respectively, both calculated for a half wave length a equal to the 
column length L. The stress cr^ refers here to buckling in the 

Z-direction (figs. 17 and 22). In particular, in the elastic range 


a l 


(a/r y ) 2 


( 122 ) 


and 



+ EC 



(124) 


where I , GI^, and EC W are the polar moment of inertia about the 

shear center S, the torsional rigidity, and the warping rigidity, 
respectively. 

In the plastic range 


a 


1 



and 


( 121 ) 


°2 = ^( EFI t + V» f|) (123) 

where F is given by equation (2) and E-fc is the tangent modulus. 

Values r Q and y 0 in the above equation for cr cr are the pr Lar radius 
of gyration about the shear center S and the distance between shear center 
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and center of gravity, respectively (figs. 21 and 22). This buckling 
stress a cr governs only if it is smaller than the column buckling 
stress for buckling in the Y-direction (figs. 17 and 22), for which 
case the interaction is practically always negligible. 

If the web of a T-section is substantially rotationally restrained 
by the flange the critical stress has to be calculated from equa- 
tions (178), (149), (119), and (189). For p = ^ (fig. l) 


a l 

a cr - 0 r i + a 2 _ 5 °2-3 


(178) 


where 


2-3 2[ 


rryyp + °3 ‘ +( a 2x + V 3 ) ' 

]f^2 + a 3 - *( a 2w + V 3)] 2 " 4(1 ' ^ )CT 3( a 2 - ^ 3 a 2 w) * 


while 


a n = 


n 2 E 


(a/ r v)‘ 


(149) 


( 122 ) 


■ 4K + 


°2V * I^( GI t» 


+ EC ™ ? 


'3 ' J + °- 12 * 2 + °' 5T ) 


^ = 


( 

1 /P 2 ) 

1 + 0.57 

( 1 /P 2 

') + 0.125P 2 + 0.57 


(124) 

(l4o) 

( 120 ) 

(136) 


^ - -'-pw/^-p 


(148) 
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The second subscript w in equations (l4o) and ( 148) refers to the web 
alone . 

In the plastic domain 


CT 1 


* 2 E t 

(a/r y ) 2 


( 121 ) 




P\ 


EFI t + E t C w 


(123) 


a 3 = + °- 12 5D3 2 + 1.84F - 0 . 25By 


(125) 


while instead of equations (l40) and (136) 


' I^( EF h» + Vw j 


q = (A/P 2 ) + 1.8UF - 0.25B 

(a/P 2 ) + 0.125DP 2 + 1.8UF - 0.25B 


All plastic values E t , A, B, D, and F refer to the actual buckling 
stress a cr . The latter four values are given in equations (2). 

The stress a cr from equation (178) has to be calculated for that 
half wave length a = pb for which it is minimum. This is practically 
the half wave length for which 0 o o is minimum. The value 

\ Lr /min 

thus obtained governs up to about a value p = a/b = Po (fig. l) . 

For p > p 2 


a cr 


a 




a l-2 


o' 2 2 


(119) 


76 


NACA TN 2640 


where 


(189) 

^ LT \° O I ~o /4 

vhile G]_, 0 2 , and cr^ are given by equations (122), (124), and (120) 

or equations (l2l), ( 123 ), and ( 125 ), for the elastic or plastic domain, 
respectively, and are calculated for a = L. The symbols r 0 and y Q 
are the polar radius of gyration about the shear center S and the dis- 
tance between shear center and center of gravity, respectively (fig. 21). 

EXPERIMENTAL INVESTIGATION 

STRESS-STRAIN TESTS 
Description of Specimens 


1-2 ofi - (v 2 /r 2fl 


c x + a 2 - /( a l - a 2 f + { ka l°^o 2 / r o C 


Square tube specimens used for this series of tests are 61S-T6 alu- 
minum alloy with the following nominal dimensions: 

2 by 2 inches by O.O 63 inch designated B 

2^ by 2^- inches by 0.047 inch designated "D" 

All tubes are special drawings of the Aluminum Company of America with 
square corners and slight thickening of the walls on the inside near 
the comers. 

Deviations from flatness, straightness, and squareness are well 
within tolerable limits, increasing as expected with wall width and 
diminishing wall thickness. The "D" tubes, in particular, show maximum 
deviation from squareness of about 5°, from flatness equal to the wall 
thickness, and a twist not exceeding 5° per 15 feet of length. One- 
third of the specimens are thus affected. Deviation from straightness 
is negligible. Variation from nominal wall thickness for all tubes is 
0.0010 to -0.0040 inch. 


Instrumentation 


In order to determine the compressive stress-strain characteristics 
of the tube material, it was necessary to prevent premature local buckling 
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of the tube wall. The walls were supported by blocking inside and out- 
side^ such that the unsupported portion of wall had a b/t ratio not 
exceeding 12.5- Ascertainable buckling was prevented at least as far 
as a strain corresponding to a secant modulus of 0.7E* 

The external blocking arrangement consisted of three square clamping 
frames which held vertical steel supporting blocks one against each face 
of the tube, as shown in figure 23 . 

Small internal clearance of the tubes made it necessary to design 
a special expanding fixture operable from the ends of the specimens. 

The device, consisting of two supporting blocks and the screw-driven 
wedge system, is shown in figures 2k and 25* The range of expansion is 
about l/4 inch, requiring the use of auxiliary blocks for the larger 
tubes. Figure 23 shows the expansor as applied to a specimen. 

Strain measurements were made with SR-4 electrical resistance 
strain gages, type A-l, in connection with standard Wheatstone bridge 
strain recorders of both Baldwin and Young manufacture. Eight gages 
were used, two to a face outside the supporting blocks, located near 
the corners at the midlength of the specimen. 


Test Procedure 

All stress-strain specimens were 8 inches long to avoid end effects 
and to provide a convenient size for handling. Nearly perfect flatness 
of ends was obtained by squaring and sanding of the sawed specimens on 
a disk sander, followed by hand-lapping on a surface plate with oil and 
emery. 

The internal expansor, slightly shorter than the specimen, was 
inserted first with necessary auxiliary blocking, and centered on the 
length of the tube. All block surfaces contacting the tube were 
lubricated with medium-weight cup grease to avoid frictional restraints. 

The external blocking was next applied. The steel supporting blocks, 
lubricated with cup grease, were centered vertically and laterally, 
supported at the base on sponge-rubber pads, and held in place by the 
center clamping frame. All blocking was then drawn up to the tube, a 
light seating load was applied, and the other two clamping' frames set 
in place. 

Preliminary tests were made to check the effect of varying clamping 
pressures. On the basis of these tests, it was decided that slight 
clamping pressure corresponding to a quarter-turn tight on all screws 
would not be detrimental. More than this showed definite gage "lag" in 
repeated load cycles, while less resulted in premature local buckling. 

The pressure so selected was used in all tests. 
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Centering the specimen, that is, providing uniform stress distribu- 
tion, was perhaps the most difficult and persistent problem encountered 
in this series of tests. To correct for nonparallelism of ends and/or 
machine heads, use was made of tissue paper shims 0.0015 inch thick 
slipped between the upper machine head and the corners of a hardened- 
steel bearing block on the upper end of the specimen. Shims were applied 
or relocated until strain readings on the eight gages showed a total 
high-to-low deviation of less than 3 percent. Usually it was possible 
to hold this to less than 1~ percent at each of three widely separated 

loads in the elastic range. 

Beyond this point no unusual problems were encountered. During the 
stress-strain tests precautions were taken to observe closely beginning 
of buckling and adequacy of the blocking and to keep a running load- 
strain curve . 


Summary of Stress -Strain Data 

Eight stress- strain specimens involving the two sizes of tube have 
been tested. Representative stress-strain curves are presented in 
figure 2 6 . Comparison of the test results showed consistent and similar 
characteristics for the two sizes of tube, with negligible deviation 
from the average for each series, a yield stress of about W -,000 psi, 

and an elastic modulus of 10.7 X 10^ psi. 

In the course of testing, 23 stress-strain tests on three other 
sizes of similar tubing were run, and they gave results completely- con- 
sistent with those of the ”B” and "D" series. 

As a further check on consistency of results, the Ramberg-Osgood 
formula (reference l8) was used for comparison with the experimental 
data and found to agree almost exactly with the average of the "B" and 
"D" series. In addition, the curves of worst deviation from the average 
were calculated and compared with the group average for stress-strain 
characteristics, tangent modulus, and secant modulus. In all cases the 
agreement was quite favorable. 


COLUMN TESTS 

Description of Specimens 


The two sizes of square tube specimens used for this program are 
described under ’’Stress -Strain Tests” and were chosen so as to have the 
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critical plate buckling stress fall in the plastic range for one series 
("B") and in the elastic range for the other ("D"). Specimens were 
selected to cover a wide range of L/r ratio, encompassing the interaction 
specimens) that is, those specimens for which the critical plate -buckling 
stress and primary column-buckling stress are identical. The interaction 
specimens in this instance are those designated I in tables 2 to 4 and 
have a ratio L/r = 45-95 for the "B" series and L/r = 85 . 80 for the 
"D" series. Including duplicate and triplicate specimens, l4 "B" columns 
and 7 "D" columns were tested, as given in tables 2 to 4. 


Procedure for Determining Effective Length of Interaction 
Column with Box Section 

(1) For series "D" specimens for which plate buckling stress 02 is 
in the elastic range: For the interaction column. 


n 2 E + 

(L/r ) 2 



from whence L/r = 1.65b/t where b is the plate width center-to- 
center of adjacent plates. 


(2) For series "B" specimens for which plate buckling stress a 2 is 
in plastic range: The length of the interaction specimen should" be such 

that 


a l 


it 2 E. 

(L/r ) 2 


a 2 


(a) Given the stress-strain curve, compute the curve for plastic 
buckling stress cr 2 p as follows. Assume several values of cf 2 p and 
for each calculate the parameters E^, E s , and 

a 2p/ a 2e = 0.^55(V^D + B + 2F), where is the elastic plate buckling 

stress. (See reference 4. table 5*) The parameters A, D, B, and F 
are given by equations ( 2 ). Further, compute for each selected value of 
cx 2 p the corresponding value of a 2e /E according to the relationship 

a 2 p cr 2 -p 

t t — . This establishes a point of the required curve 
E ( a 2 p/ cr 2 e) E 

<? 2 p against CT 2 e/ E f° r each value of crgp 


selected. 
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(b) Having found the curve for cr 2p , calculate 



With this value as abscissa, pick the corresponding cr 2p from the curve, 
and the corresponding E^ to the stress-strain curve at that stress. 

Then calculate the required L/r ratio as defined above. 

(c) The curve for c 2p ma y be plotted on the same coordinate system 
as the stress-strain curve, as shown in figure 27. 

It will be noted that E^ and a 2p are interdependent, which 
necessitates the trial-and-error procedure as given in order to find the 
values satisfying the interaction criterion. 

Since c 2e /E is immediately known and the curve for ff 2p is 

usually near to the stress-strain curve, it is necessary to select only 
the values of "that narrowly encompass the stress corresponding 

to this value of a 2e/^’ 

(d) The specimen length may be determined by correcting the effective 
length just computed for the particular end-support conditions used. 


Instrumentation 

In order to study properly the interaction of local and Euler 
buckling, it was necessary to develop means of measuring separately the 
two types of buckling. It was felt that electrical resistance strain 
gages on the column faces would not accomplish this by virtue of 
difficulty in "sorting out" the proportionate effects of the two types 
of deformation. In terms of mechanical gages, it was recognized that 
a device to measure plate deflections due to local buckling must be 
attached to and "ride" with the column during primary deflection in 
order to exclude the effects of the latter. The primary column deflec- 
tions are measurable by any one of several simple devices referred to 
the ends of the column. 

A collar to fit on the column and carry a local buckling gage without 
affecting the local buckling characteristics has been developed. 

The local buckling gage developed for these tests shown in figure 28 
is based in principle on the type of gage described in reference 20 and 
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has the advantage of being independent of buckle location. It makes use 
of a suspended blade in contact with the column face and measures the 
blade movement resulting from buckle formation by means of dial gages. 
(The gages were Federal No. C21 with 0.0001- in. divisions and 0.3- in. 
range and were full- jeweled. ) The blade suspension is 1/2 by 0.002-inch 
stainless spring steel, combining high resistance to breakage and 
negligible blade restraint. 

Gages of this type were applied (at two different levels} on the 
column faces parallel to direction of primary buckling in order to avoid 
effect of primary curvature of the column. 

Euler deflections were measured by means of 0.001-inch dial gages at 
midlength and both quarter points, in combination with a 0.01- inch 
division scale at midlength read simultaneously for correlation with 
reset readings of the dial gages and for large deflections beyond the 
range of the dial gages. 

Columns were supported at the ends by knife edges, with appropriate 
adjustment for centering. Carboloy knife edges, with corresponding flat 
bearing surfaces, were loaned by the National Bureau of Standards. 

A general view of the testing arrangement is given in figure 29. 


Test Procedure 

The ends of all columns were squared, sanded, and hand-lapped in 
a manner similar to that for the stress-strain specimens. 

The centering procedure made use of both strain readings and lateral- 
deflection readings. Eight Tuckerman optical strain gages were applied 
to the corners of the column at two stations, generally 12 to 15 inches 
outside the upper and lower quarter points, and always at a distance 
from the ends of at least twice the plate width to avoid end effects. 

The two sets of strain readings made a simple matter of determining 
required centering adjustments and in particular of deciding whether 
eccentricities of the column ends were in the same or opposite directions. 

Final centering, at about two-thirds of the predicted critical load, 
was done by adjusting the column ends until lateral deflections at mid- 
length and both quarter-points were negligible, and then making a final 
check of strain distribution. Centering by deflection proved to be 
considerably more accurate than strain readings for the final adjustment. 
Differences in the average strain on opposite faces of the column were 
held to a maximum of 1 percent in the direction of primary buckling and 
two percent in the other less critical direction. 
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Local buckling gages were then applied to the centered column, and 
the test proceeded by reading local and Euler deflections at each load 
increment to failure . 


Evaluation of Experimental Data 

Ultimate stress for each column is given in tables 3 and k and also 
is plotted in relation to the tangent-modulus column curves in figures 30 
and 31 * 

The method of obtaining the tangent modulus from the stress-strain 
curve may be of interest here. A semitransparent mirror with a one-to- 
one transmission-reflection ratio was used. It was placed perpendicular 
to the curve at the desired point and perpendicularity was checked by 
casting the reflected image of a short portion of curve in front of the 
mirror onto a short portion of the curve seen through the mirror. A 
line - the normal to the curve - was then scribed along the face of the 
mirror. The method proved to be more accurate than drawing the tangent 
by eye, especially in the sense that several operators with little 
experience could obtain precisely the same value of modulus. 

Local buckling data were evaluated by means of the "top-of-the- 
knee" method as developed by the NACA (reference 19) and critical plate- 
buckling stresses thus obtained are plotted on the "D" -series column 
curve, figure 31 • Typical curves of load against local deflection for 
the "D"-series columns are shown in figure 32 . No local buckling data 
were obtainable for the "B" series, since the critical plate stress was 
in the plastic range with the result that beginning of plate buckling 
was immediately followed by complete collapse of the column. Consequently, 
the ultimate load for the column coincided with the plate buckling stress, 
and no postbuckling strength was indicated. This observation is in 
accordance with results given in reference ( 21 ). In contrast, the 
"D"-series columns had considerable postbuckling strength as shown in 
table 4 and figure 31* 


CONCLUDING REMARKS 


The following results were observed from an investigation made to 
determine the interaction of column and local buckling in compression 
members . 

1. The test results do not show any noticeable interaction effect 
Por the square tubes tested. This seems to substantiate the theory, 
which shows the order of stress reduction due to interaction effect to 
be so small that it cannot be differentiated experimentally from the 
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reduction of concentric buckling stress occurring as a result of inherent 
eccentricities. 

2. The interaction effect is negligible for box sections, as 
indicated by both theory and experiment. As stated in the paper, the 
same conclusion applies to the common sizes of I-, H-, and channel 
sections, but not to sections for which torsional instability is an 
important factor, such as T and angle shapes. 

3- Tests confirm the known fact that appreciable postbuckling 
strength beyond the critical plate stress is possible when the latter 
is well within the elastic range. Conversely, when the critical plate 
stress is in the plastic range, complete collapse of the column accom- 
panies beginning of plate buckling. 


Cornell University 

Ithaca, N. Y. , March 23, 1951 
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TABLE 1.- VALUES OF a FROM THEORETICAL INVESTIGATION 


CD 


a 

( cm) 


0 2 

03 

ct 2 v 

°l-2 

ct 2-3 

a cr 

from 

equation 

(33) 

T 

°CT 

from 

equation 

(177) 

a cr 

from 

equation 

(178) 

a cr 

from 

equation 

(208) 

a cr 

from 

equation 

(119) 

Unit 

50 







2,440 

34,650 





kg/cm 2 

psi 

58.33 

304,670 

4,325,000 

6,030 

85,600 

2,553 

36,220 

1,323 

18,800 


2,410 

34,240 

2,375 

33,750 

2,405 

34,170 

2,391 

33,950 



kg /cm 2 
psi 

62 

269,800 

3,830,000 

5,917 

84,000 

2,607 

37,000 

1,267 

18,000 

5,889 

83,600 

2,427 

34,460 


2,421 

34,380 

2,405 

34,150 



kg/ cm 2 
psi 

87.5 

135,580 

1,925,000 

5,483 

77,900 

3,333 

47,350 

1,042 

14,800 

5,444 

77,250 

2,755 

59,120 

2,755 

39,120 

2,741 

38,920 

2,700 

38,350 



kg/cm 2 

psi 

135 

56,900 

807,500 

5,230 

74,300 

5,907 

84,000 

910 

12,910 


3,533 

50,200 


3,477 

49,400 

3,326 

47,220 



kg/cm 2 

psi 

175 

33,870 

480,400 

5,156 

73,300 

9,043 

128,300 


4,961 

70,500 


3,855 

54,750 





kg/cm 2 

psi 

350 

8,467 

120,100 

5,067 

72,000 

32,600 

463,000 

830 

11,790 

4,289 

60,900 

4,613 

65,500 

4,030 

57,250 



4,019 

57,100 

3,790 

53,850 

kg/cm 2 

psi 

620 

2,699 

38,300 

5,052 

71,800 

99,800 

1,417,200 

819 

11,630 

2,341 

33,250 





2,325 

33,000 

2,287 

32,500 

kg/cm 2 

psi 

700 

2,117 

30,050 

5,050 

71,800 

127,000 

1,804,000 

817 

11,600 

1,911 

27,170 

4,925 

70,000 

1,880 

26,700 



1,904 

27,040 

1,883 

26,750 

kg/cm 2 

psi 
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TABLE 2.- PROPERTIES OF TUBE SPECIMENS 

^Series "B" tubes are 2 by 2 by 0.062 in. with av. r = 0.792 
in.; series "D" tubes are 2^ by 2^- by 0.047 in. with 
av. r = 1.001 in.l 


Tube 

Col- 

umn 

Cut 

length 

(in.) 

Weight 
( grams ) 

Area 
(sq in.) 

Corrected 

free 

length, 

L 

(in. ) 

L/r 




Series ' 

’B" tubes 


"B"-l 

SS -1 

8.05 

171.7 

0.481 

11.82 

14.95 

"B"-2 

SS-2 

8.11 

169 .I 

.470 

11.88 

15.00 

"B "-1 

ES-1 

11.38 

238.9 

.474 

15.23 

19.25 

"B"-2 

ES-2 

11.39 

238.6 

.473 

15.24 

19.25 

"B"-l 

S-l 

17.48 

368.7 

.476 

21.45 

27.15 

"B"-2 

S-2 

17.49 

366.6 

.474 

21.46 

27.15 

"B"-l 

M-l 

26.00 

551.8 

.478 

30.02 

37.95 

"B"-l 

M-2 

26.00 

549.2 

.476 

30.02 

37.95 

"B"-2 

1-1 

32.34 

677.3 

.473 

36.37 

45.95 

"B"-2 

1-2 

32.36 

677.2 

.473 

36.39 

45.95 

"B"-2 

1-3 

32.36 

676.2 

.472 

36.39 

45.95 

"B"-3 

ML-1 

36.35 

771.4 

.479 

40.38 

51.00 

"B"-3 

L-l 

40.23 

853.0 

.479 

44.26 

56.00 

"B"-3 

L-2 

40.23 

850.2 

.477 

44.26 

56.00 




Series " 

D" tubes 



"D"-4 

SS-1 

17.31 

346.2 

0.452 

21.25 

21.23 

"D"-l 

S-l 

38.88 

780.0 

.452 

43.03 

43.00 

"D"-2 

RM-1 

67 .OO 

1349.7 

.455 

71.09 

71.02 

"D"-4 

MS-1 

72.16 

1446.0 

.452 

76.25 

76.17 

"D"-3 

1-2 

81.78 

1635.2 

.452 

85.87 

85.78 

"D"-l 

1-1 

8l.8l 

l64l.O 

.453 

85.90 

85.81 

"D"-5 

L-l 

123.49 

2545.0 

.465 

127.58 

127.45 
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TABLE 3.- TEST RECORD OF SERIES "B" TUBES 


Dimensions : 


2 by 2 by 0.062 in. 


Designation 

Theoretical 

Experimental 

P ult/ P cr 

Tube 

Column 

Plate 

a cr 

Column 

a cr 

P cr 

P ult 

a ult 

"B"-l 

SS-l 

37,250 

43,500 

17,900 

17,850 

37,150 

0.997 

"B”-2 

SS-2 

37,250 

43,500 

17,500 

17,700 

37,650 

1.010 

"B"-l 

ES-1 

37,250 

42,700 

17,630 

17,900 

37,800 

1.015 

m B m -2 

ES-2 

37,250 

42,700 

17,600 

17,400 

36,800 

.988 

”B M -1 

S-l 

37,250 

40,900 

17,720 

17,750 

37,270 

1.001 

"B"-2 

S-2 

37,250 

40,900 

17,630 

17,450 

36,810 

.988 

M B M -1 

M-l 

37,250 

38,500 

17,800 

17,500 

36,600 

.983 

"B”-l 

M-2 

37,250 

38,500 

17,720 

17,250 

36,200 

.972 

M B M -2 

1-1 

37,250 

37,250 

17,600 

16,900 

35,750 

.960 

M B M -2 

1-2 

37,250 

37,250 

17,600 

17,400 

36,800 

.988 

M B M -2 

1-3 

37,250 

37,250 

17,590 

17, 150 

36,350 

.975 

"B"-3 

ML-1 

37,250 

36,500 

17,470 

17,150 

35,800 

.980 

m B m -3 

L-l 

37,250 

33,600 

16,100 

15,400 

32,200 

.959 

m b m -3 

L-2 

37,250 

33,600 

16,030 

15,400 

32,350 

.963 


TABLE 4.- TEST RECORD OF SERIES "D" TUBES 
dimensions: 2 ^ by 2 -i by 0.047 infj 


Desig] 

aation 

Theoretical 

Experimental 

Ratio of experimental 
to theoretical values 

Tube 

Column 

Plate 

a cr 

Column 

CT cr 

P cr 

P ult 

a ult 

°cr 

(a) 

■^ult/^cr 

Plate 

a cr/ a cr 

"D"-4 

"D"-l 

"D"-2 

"D"-4 

"D"-3 

"D"-l 

"D"-5 

SS-l 

S-l 

RM -1 

MS-1 

1-2 

h I-l 

L-l 

1^,325 

1^,325 

1^,325 

14,325 

1^,325 

1^,325 

1^,325 

42.000 
37,650 
20,700 

18.000 
1^,325 
1^,325 

6,500 

6470 

6470 

6520 

6470 

6480 

6500 

3020 

9670 

9200 

6800 

6525 

6230 

6050 

2900 

21,390 

20,350 

14,920 

14,420 

13,790 

13,330 

6,240 

14,290 

14,200 

13,850 
l4 , 260 

1.491 

1.420 

1.042 

1.007 

.963 

.930 

.960 

0.997 

.991 

.967 

.995 


Experimental plate a cr 
^Two tests on same column. 


as determined by NACA top-of-knee method. 
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Figure 1.- Diagram of "buckling stresses plotted against ratio of half 
wave length to web width 3 = a/b. 
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(a) Cross-sectional sketch of column. 



(b) Cross-sectional sketch of web showing interaction of column and 

plate buckling. 



(c) Sketch of column showing axis notation. 


Figure 2.- I-section column showing various notations. 


Buckling stress, cr cr , 


92 


NACA TN 2640 



Figure 3.- Calculated buckling stresses of T-section of figure 16 for 

several half wave lengths a. 
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Figure 
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t l’ 


v" 


£ y j 

>. 
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(a) Cross-sectional sketch of column. 

y 



(b) Cross-sectional sketch of web shoving instability. 

4.- Asymmetric column supported at both unloaded edges by flanges 
of different width and thickness. 



Figure 5.- Sketch of lover half of flange. 
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Figure 6.- Sketch of T-stiffener and column made of T-sections. 



(a) Sketch of column showing various notation. 






W 1 




(h) Buckling as column without distortion of cross section. 


V 


\ 




f^2 
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(c) Buckling as plate with distortion of cross section. 
Figure 7.- Sketch of I-section with narrow and thin flanges. 
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Figure 8.- Curves of against 0 in elastic range for I-section 
columns with flanges of 4- and 2-inch width. 
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Figure 9-- Diagram of square tube or box section. 





Figure 10.- Deflective forces 


-ta ^ 
dx 2 
section. 


caused by deflection of box 
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Figure 11.- Square tube with plates buckling symetrically with respect 

to vertical and horizontal axes. 



Figure 12.- Square tube showing alternate mode of web buckling. 


rolcr 1 
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Figure 13.- H-section where width of flanges is such that they are 
rotationally restrained from buckling by web. 



t 



b 



Figure 14.- Sketch of T-section showing notations. 
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Figure 15.- T-section with small rotational restraint of web. 



Figure 16.- T-section where flange has relatively large torsional 
rigidity and is much more stable than the web, so that web is 
substantially rotationally restrained. 
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Figure 17.- Deflection of cross section in figure l 6 consisting of 

translation w-^ rotation with respect to shear center S, and plate 
buckling of web with maximum deflection W 3 . 



Figure 18.- Deflections of case (2) (w 2 ) and case (3) (w^). 
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Figure 19*- Deflections of case (l) and of combined case (2-3). 


z 



Figure 20.- Distortion of cantilever t>eam under uniform load q. 
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Figure 21.- Flexural and torsional buckling of T-section showing 
displacement CC ' of center of gravity C. 



y 


Figure 22.- Flexural and torsional buckling of angle section. 
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Figure 24.- Internal expansor for tube specimens, assembled 
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Figure 25.- Internal expansor showing individual parts. 
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Figure 26.- Typical stress-strain curves for square-tube specimens. 
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Figure 27 • - Curve for 


CT 2p 


plotted on same coordinate system as stress- 
strain curve. 
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Figure 28.- Local buckling gage and supporting collar. 
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Figure 29.- General view of test arrangement 
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Figure 31- - Test results for series "D" specimens. 
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